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ABSTRACT 

The main goal of this paper is to compare the relative importance of destruction by tides, vs. 
destruction by mergers, in order to assess if tidal destruction of dwarf galaxies in clusters is a viable 
scenario for explaining the origin of intracluster stars. We have designed a simple algorithm for 
simulating the evolution of isolated clusters. The distribution of galaxies in the cluster is evolved using 
a direct gravitational N-body algorithm combined with a subgrid treatment of physical processes such 
as mergers, tidal disruption, and galaxy harassment. Using this algorithm, we have performed a total 
of 227 simulations. Our main results are (1) destruction of dwarf galaxies by mergers dominates over 
destruction by tides, and (2) the destruction of dwarf galaxies by tides is sufficient to explain the 
observed intracluster light in clusters. 

Subject headings: cosmology — galaxies: clusters — galaxies: dwarfs — galaxies: interactions — 
methods: numerical 
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1. INTRODUCTION 

1.1. Dwarf Galaxies 

Dwarf galaxies (DGs) are the most numerous galaxies 
occurring in the Universe. A ma j ority of galaxies in 
the local group are DGs (jMateal L1998'). Also DGs 
comprise 85% of the Local Volume galaxy population 
{D < 10 Mpc, iKarachentsev et all |20M) . and have 



Coma (|Thompson & Greeorvl 


19931: Bernstein et al. 


19951), Virgo ISandage et al. 


1985': 'ImDcv et al. 


19881: iPhilliDDS et al. 1998: Lee et al. 20031. Fornax 


(Bothun et al.l 11991 


I iDrinkwater et al.l 120031. Centau- 


rus (iMieske et al.l 


200711, and several galaxy groups 



(iKarachentseva et al.l 



199H ICotc ct all fl9M 



Cellone fc B uzzoni l2005h . 
density ~ 40 times that 



ICarrasco etHI 120011 : 

DGs may have a space 
of bright galaxies in the 



Universe (jStavelev-Smith. Davies. fc Kinmanlll992l l. 

DGs are defined as low-mass (10^ - 10^ Mq) galaxies 
having an absolute magnitude fainter than Mb ^ — 16 
mag, or My ^ — 18 mag ( Grcbcl 2001), have low sur- 
face brightness and low metallicity. Their small stellar 
fraction and very low luminosities make them the hard- 
est galaxies to detect. They are believed to be the sin- 
gle systems with the largest proportion of dark-matter, 
and have a correspondingl y high ratio of dark to lumi- 
nous mass (e.g.. lC6te. Carignan. fc Freemanll2000( l. with 
M/L ratios as high as that of galaxy groups and poor 
clusters. Due to their smaller masses and gravitational 
potentials, DGs are less able to retain their gas as com- 
pared to more massive galaxies, and in a clustered en- 
vironment, the DGs are more likely to be disrupted by 
galactic encounters and environmental effects. 

In the hierarchical clustering scenario of structure 
formation in the Universe (e.g.,JWhite &: FrenkI 119911 : 
iKauffmann. White, fc Guiderdonill993f l. the accretion of 
DGs causes the build up and growth of massive galax- 
ies and large-scale structures. There exists a deficiency 

^ Departement de physique, de genie physique et d'optique, Uni- 
versite Laval, Quebec, QC, Canada 



in the number of observed low-luminosity DGs (dis- 
crepancy more than one order of magnitude) as com- 
pared to the large number of th eoretically predicted low- 
mass dark mat ter halos (e.g ., 'Trentham fc Tullvl l200l : 
iTrcntham, Tu Uv. fc Mahdavi 2006). This is recognized 
as a problem for cold dark matter theory, and the likely 
solution to this problem involves energy feedback from 
stellar evolution. Dense clusters are observed to con- 
tain a larger number of low-luminosity DGs per high- 
luminosity giant galaxy wh en compared to the field. 
iTrentham fc HodgkinI ()2002l ) found that the Virgo clus- 
ter contains ~ 2.5 times more dwarf galaxies per giant 
galaxy when compared to the Ursa Major cluster. These 
imply that dwarfs are more common relative to giants in 
dense environments than diffuse ones. 

The dynamical evolution of galaxies in a cluster is in- 
fluenced by several mechanisms. There are two types 
of tidal interactions: t he tidal force s due to other (mas- 
sive) cluster galaxies (|Gnedinll2003| l. and the tidal field 
resulting from the overa ll cluster potential (Merritt 198^ 
iBvrd fc Valtonenlll990l l. There can be collisions between 
the galaxies themselves due to their motion, sometimes 
resulting in mergers. Of particular interest is the oc- 
currence of mul tiple hig h-velocity encounters between 
cluster galaxies (jRichsto nc 1976), a phen ome non which 
is termed "galaxy harassment " (iMoore e t al. 1996). A 
phenomenon like tidal stirring ijMaver et al.y2001,. where 
tidal shocks strip DGs) has a more pronounced effect on 
the less massive galaxies. There can also be stripping 
within galaxy group s and protoclusters accreting onto a 
cluster (|Mihosl l2004l. 

Works are found in the literature on the interaction 
of cluster environment with DGs in the cluster. Study- 
ing the core of the Fornax cluster, iHilker et al.l (119991 1 
examined a scenario in which DGs are accreted and dis- 
solved at the cluster center. They found that the infall 
of DGs can largely explain many properties, but there 
are probably other physical processes occurring simulta- 
neously, like the evolution of the more mass ive galaxies 
by stripping and merging. iMori fc Burkei^ ([2000) per- 
formed analytical and hydrodynamic studies of the in- 
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teraction between the ICM of a cluster and an extended 
gas component of DGs surrounded by a cold dark mat- 
ter halo. They found that the DG halos lose their diffuse 
gas rapidly by ram-pressure stripping in a typical cluster 
environment. 

1.2. Intracluster Light 

The diffuse intracluster light (ICL) observed in clus- 
ters of galaxies is produced by stars, usually of low 
surface brightness, located outside individual galaxies 
but within the cluster and associated with the clus- 
ter potentia l . Th e first mention of IC light was made 
by IZwickvl (jl951f ). who observed extended irregular 
regions of stars and low surface-brightness matter in 
the intergalactic spaces of the Coma cluster. Sev- 
eral observations hav e since detected diffuse light in 
galax y clusters (e.g., IVilchez-Gomed Il999t lArnaboldH 
|2004 for reviews). Diffuse ICL has been observed in 
many galaxy cluster systems (^e.g.. iGonzalez et 311120051 : 



Krick et a l. 2006), including non-cD galaxy clusters (e.g. 
Feldmeie i_ et al. 200 4a'). Observations of diffuse ICL in 
Virgo (|Mihos et al]|2005f ) reveal complex substructures 



in intricate web patterns including long tidal streamers, 
small tidal tails, and intergalactic b ridges. The idea o f 
IC globular clusters was proposed by IWest et all (|1995D . 
Later on, distinct IC stars were observed, including glob- 
ular clusters, red giant stars, and SN la (|Gal-Yam et al.l 
|2003|). In dividual IC stars have been fou nd in Virgo 



Ferguson e t al. 1998; Arnaboldi et al.ll2003l ). and Coma 
Gerhard et al.i i2005j). Deep broadband imaging has 



helped to observe ICL in several other local galaxy clus- 
ters, e.g., Cantaurus, Fornax, Abell, and clusters up to 
intermediate redshifts. The origin and evolution of the 
IC stars and diffuse light are not well constrained at 
present. 

Arguably these field stars are gravitationally stripped 
from their parent galaxies when the galaxies in a clus- 
ter interact dynamically with other galaxies or with the 
cluster potential. Continuous accumulation of ICL is 
now believed to be an ubiquitous feature of evolved 
galaxy clusters, with the unbound (i.e. intracluster) 
star fraction slowly increasing with time. From ob- 
servations and cosmological simulations, at z = at 
least 10 — 20% of all st ars in a cluster are unbound 
to any one galaxy (e.g., lAguerri et al.l l2005f ). As an 
upper limit, ICL constitute « 10 — 40% of a cluster's 
luminosity. In observations and simulations, the frac- 
tion of stars in ICL increases with mass of the clusters, 
and increases with density of environment: from loose 
groups (< 2%, Castro- Rodriguez et al. '2003"), to Virgo- 
like (10%, Fcldmcicr et al. 2004b; Zibctti ct al. 2005) and 
rich clusters (^ 20% or higher. ITvson fc Fischer! 119951 : 
iFeldmeier eralll2004al: iKrick fc Bernsteinll2007l l In the 
cores of dense and rich clusters ( hke Coma) the local IC L 
fraction can be as high as 50% ([Bernstein et al.lll995( l. 

The most popular formation mechanism of IC popula- 
tion is stripping of stars from cluster galaxies by grav- 
itational tides, fast encounters between galaxies, and 
tidal interact i ons between colliding and merg i ng galaxies 
(piiie? 1983'; 'Richston e fc MalumuthI [Ml IWeil et all 
[1997; Gregg fc West 199^. Another contribution to the 
IC population can come from the infall of galaxy pro- 
toclusters containing stars which are al ready unbound. 
Studies of IC stars in the Virgo cluster ([Feldmeier et al.l 



l2004bf) imply that bulk of the IC stars come from late- 
type galaxies and dwarf galaxies, and that the IC stars 
form by tidal stripping. Optic al studies of IC L in in- 
termediate redshift clusters (Zib etti et al.|[2005[ ) support 
that ICL is produced by stripping and disruption of 
galaxies when they pass through cluster centers. A dom- 
inant mechanism of ICL formation is believed to be tidal 
stripping during the hierarchical assembly of clusters. 

Several authors have studied the origin of the dif- 
fuse ICL in clusters of galax ies using numerical simula- 
tion. iNapolitano et al.1 (|2003f ) performed high-resolution 
N-body ACDM simulations of a Virgo-like cluster to 
study the velocity and clustering properties of the IC 
stellar component at z = 0. The member galaxies in 
their simulated cluster undergo tidal interactions among 
themselves and with the cluster potential, to produce 
diffuse stars in the ICM . Their results substantially 
agree with the observed clustering properties of the dif- 
fuse IC stars in Virgo. Using hydrodyna mical simulation 
of a ( 192 Mpc)^ cosmological box, iMurante et al.l 
(|2004D studied the statistical properties of IC stars in 
galaxy clusters and the dependence of the ICL prop- 
erties on cluster mass and temperature. Their simula- 
tions reveal substantial (at least ^ 10%) diffuse stellar 
light in a cluster, hidden from observations because of 
its very low surface brightness. Using high-resolution N- 
body -I- SPH simulations of a Coma-like c l uster formed 
in a cosmological context, IWillman et al.l (|2004( ) stud- 
ied the formation, evolution and properties of IC stars 
in a rich cluster. Their results indicate that both mas- 
sive and smaller galaxies contribute to ICL formation, 
with the stars stripped preferentially from the outer, 
lower-metallicity parts of a galaxy's stellar distribution. 
The resulting fraction (^ 20% at z = 0), and distri- 
bution of IC stars are in good a greement with the ob- 
served ICL properties. iSommer-Larscn et al. (2005) per- 
formed cosmological TREESPH simulations of the for- 
mation and evolution of galaxy groups and clusters, in 
order to discuss the ICL properties. They simulated a 
Virgo-like and a Coma-like cluster, and 4 galaxy groups, 
and predict several properties of IC stellar populations. 
iRudick et al.l ([20061 ) have used A^-body simulations (done 
at both low and high resolutions) of the dynamical evolu- 
tion of galaxy clusters to study the formation and evolu- 
tion of the diffuse ICL component in 3 simulated galaxy 
clusters. They found that the ICL fraction of cluster lu- 
minosity increases as clusters evolve, reaching ^ 10—15% 
in evolved clusters. 

In these numerical studies, there is always a trade- 
off between having goo d res olution o r good statis- 
tics. Naoolitano et al. (I2003D : iWillman et al.l (12004 



ISommer-Larsen et al.l ([2005D . and iRudick et al.l ([2006 
simulate either one cluster or a few clusters, so even 
though their simulations have high resolution, they have 
poor statistics, in the sense that the cluster(s) they are 
simulating might not be representative of the whole clus- 
ter po pulation. At the other extreme, IMurante et al.l 
([2004D simulate a very large cosmological volume, con- 
taining a statistically significant sample of clusters. Such 
large simulations however cannot resolve the scale of 
dwarf galaxies. Our goal is to have it both ways: achiev- 
ing good statistics while resolving the processes respon- 
sible for destroying dwarf galaxies. This is achieved by 
combining large-scale cosmological simulations with a 
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semi-analytical treatment of mergers and tidal disrup- 
tion. 

1.3. Objectives 

The main objective of the present work is to deter- 
mine if dwarf galaxies in clusters are more prone to de- 
struction by tides or to destruction by mergers. This 
determination is then used to predict the contribution of 
dwarf galaxies to the origin of intracluster stars in differ- 
ent types of cluster environments. The DGs in a cluster 
can be tidally disrupted (by the field of a more massive 
galaxy or by the background halo,) or the DGs can be 
destroyed when they merge with another galaxy. The 
impact of these two destruction mechanisms on the ICL 
is radically different. In the case of tidal disruption, the 
process contributes to IC stars in the cluster. In the case 
of merger, the DG is absorbed by a more massive galaxy, 
and there is essentially no contribution to the IC stars. 

For the first part of this project, presented in this pa- 
per, we perform numerical simulations of isolated clusters 
of galaxies, in order to examine which method of dwarf 
galaxy destruction is dominant, and how the process de- 
pends on environmental factors. We identify six possi- 
ble outcomes for our simulated galaxies: (1) the galaxy 
merges with another galaxy, (2) the galaxy is destroyed 
by the tidal field of a larger galaxy but the fragments ac- 
crete onto that larger galaxy, (3) the galaxy is destroyed 
by tides and the fragments are dispersed in the intra- 
cluster medium (ICM), contributing to the intracluster 
light, (4) the galaxy is destroyed by the tidal field of the 
background halo, (5) the galaxy survives to the present, 
and (6) the galaxy is ejected from the cluster. 

We designed a simple algorithm to follow the evolu- 
tion of galaxies in an isolated cluster. The gravitational 
interaction between galaxies is calculated by a direct N- 
body algorithm. The other physical mechanisms gov- 
erning the possible outcomes (mergers, tidal disruption, 
accretion etc.) of the simulated galaxies are treated as 
"subgrid physics," and are incorporated in the algorithm 
using a semi-analytical method. In the present work, we 
use this algorithm to simulate the evolution of isolated 
galaxy clusters, i.e. we assume that the cluster has al- 
ready formed with its constituent galaxies in place, and it 
is neither accreting nor merging. Cluster accretion possi- 
bly has a finite impact on the evolution of dwarf galaxies, 
and on the origin of intr acluster stars. In a forthcoming 
paper (|Brito et al.l 120081 ) we will address these issues by 
doing actual cosmological simulations over a statistically 
significant volume of the universe. 

The remainder of this paper is organized as follows: In 
^ we outline the numerical model for our galaxy clus- 
ters. The methodology of our simulations is described in 
SJJI The results are presented in 21 s-^d ^JSJ We discuss 
the implications of our main goals in ^ and give our 
conclusions in fJT] 

2. THE NUMERICAL METHOD 

2.1. The Basic PP Algorithm 

We treat the system as an isolated cluster consisting 
of N galaxies of mass m^, radius Si, and internal energy 
Ui, orbiting inside a background halo of uncollapsed dark 
matter and gas. We assume that the halo is spherically 
symmetric, and its radial density profile Phaio('') does not 



evolve with time (hence, we are neglecting infall motion 
that would result from cooling flows). Furthermore, we 
assume that the halo is stationary: it does not respond 
to the forces exerted on it by the galaxies, and therefore 
its center remains fixed at a point that we take to be the 
origin. 

We represent each galaxy by one single particle of mass 
rrii. The "radius" Si of the galaxy and its "internal en- 
ergy" Ui are internal variables that only enter in the 
treatment of the subgrid physics described in ii2. 41 below. 
Our motivation for using this approach is the following: 
To simulate the destruction of dwarf galaxies by tides, 
it would seem more appropriate to simulate each galaxy 
using many particles. Supposing, however, that it takes 
at least 100 particles to properly resolve a dwarf galaxy 
experiencing tidal destruction, as the galaxies in our sim- 
ulations cover 3 orders of magnitude in mass, the most 
massive ones would be represented by 100,000 particles. 
Even though the dwarf galaxies are much more numer- 
ous than the massive ones, the total number of particles 
would be above 1 million. This raises the following issues: 

• With the use of tree codes, an iV = 10^-particle 
simulation is not considered prohibitive anymore. 
However (1) our model has several free parameters, 
so we have a full parameter-space to study, and (2) 
one single cluster is not statistically significant, so 
for each combination of parameters we need to per- 
form several simulations. For this paper, we per- 
formed 227 simulations. Doing 227 million-particle 
simulations would start to be computationally ex- 
pensive. 

• We could use unequal-mass particles, so that the 
most massive galaxies would not be represented by 
large numbers of particles. This is usually not a 
good idea. N-body simulations with particles hav- 
ing wildly different masses are known to suffer from 
all sorts of instability problems, which often require 
special algorithms to deal with. The approach we 
are considering here is more practical. 

• In this paper, we co nsider isolated clu sters. In 
a forthcoming paper (jBrito et al.l l2008f ). we will 
present simulations of a cosmological volume con- 
taining at least 100 clusters. The number of parti- 
cles would then reach 100 million, and we would 
still need to explore the parameter space. This 
would be very computationally expensive. We will 
solve this problem using single-particle galaxies 
combined with a treatment of subgrid physics. The 
simulations presented in this paper can be seen as 
a test-bed for this approach. 

The relatively small number of particles in our simula- 
tions (typically less than 1,000) enables us to use a direct, 
Particle-Particle (PP) algorithm, which is the simplest 
of all N-body algorithms. We took a standard PP code, 
which evolves a system of N gravitationally interacting 
particles using a second-order Runge-Kutta algorithm. 
We modified the original algorithm to include the inter- 
action with the background halo, and we added several 
modules to deal with the subgrid physics. In this mod- 
ified algorithm, the number of particles N can vary, as 
they merge, are destroyed by tides, or escape the cluster. 
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2.2. Gravitational Interactions 
The acceleration of particle i (or galaxy i) is given by 



-G 



(1) 



where and are the positions of particles i and j, re- 
spectively, rrij is the mass of particle j, Mhaio(''i) is the 
mass of the background halo inside r = , G is the grav- 
itational constant, and e is the softening length. This 
assumes that the background cluster halo is spherically 
symmetric and centered at the origin. In our PP algo- 
rithm, this expression is evaluated directly, by summing 
over all particles j ^ i- The softening length e is cho- 
sen to be smaller than the initial radius of the smallest 
galaxies (see §3.2 below for the determination of the ini- 
tial radius). 

We evolve the system forward in time using a second- 
order Runge Kutta algorithm. The timestep At is calcu- 
lated using 



At = min (Ai); , (At)^ = min 



e 



1/2- 



(2) 

where is the velocity of particle i, and we take the 
smallest value of {At)i to be the timestep At. 

2.3. The Cluster Halo Density Profile 

We consider two different types of density profile of the 
background halo of a cluster, phaio(f): the /3 profile, and 
the NFW profile. 

In the first case, we assume that the dark matter 
in the background halo follows a similar density dis- 
tribution as the observed intracluster gas. A single 
/?-model (isothermal ) density profile is used for the 
gas (e.g. iKind [19621: iCavaliere fc Fusco-Femianol [T976I : 
iMakino. Sasaki, fc Sutolll998f ). 



Pgas (r) = Po 1 + (r/rc) 



-3/3/2 



(3) 



where, po is the central density, and Tc is the core 
radius. The values of Pn , Tc and (3 are taken from 
iPiffaretti fc Kaastral (|2006D . which gives the gas den- 
sity parameters for 16 nearby clusters. The halo den- 
sity is then obtained by scaling the gas density with the 
universal ratio of matter (dark -I- baryonic) to baryons, 
Phaio = Pdm + Pgas = Pgas^M/^b, whcrc and fit 
are the present matter (baryons -I- dark matter) density 
parameter and baryon density parameter, respectively. 
This assumes that the cluster baryon mass fraction fol- 
lows the cosmic val ue of r^b/r^Af which is expected to be 
genera l ly tru e (e.g.. IWhite et all 1 19931 : 1 Allen et al.ll2002l : 
lEttoril l2003l ). although precise estimations of cluster 
baryor i content have shown deviations from the universal 
value ijGonzalez et al.ir2007l . and references therein). 

In a second case, we consider that the distribution of 
gas and dark matter in the background halo both follow 
analytical models of the dark matter density having a 
functional form 



(4) 



(r/rs) (1 -I- r/rs) 

avarro. Frenk. fc Whit3fT99l . Here, p^ is a scaling 
density, and is a scale length. The NFW profile is 



often parametrized in terms of a concentration parameter 
c. The parameters ps and rs are then given by 



200c3perit(^)/3 



25i?2(z)c3/7rG 



ln{l + c)~ c/ (l+c) ln(l + c)-c/(l+c) 



' 200 



, (5) 



(6) 



max] 
halo J 



where Pcrit(-z) ~ 3H'^(z)/8ttG is the critical density at 
formation redshift z, and r2oo, the virial radius, is the 
radius of a sphere whose mean density is 200pciit (200 
times the critical density of the Universe at the epoch 
of formation). After scaling, the halo density profile is 

Phalo = PdM^m/ {^M - ilfc). ^ 

Once we have chosen a particular density profile, the 
density is integrated to get the background cluster halo 
mass as ^ 

Mhaio(r') = / ■iTrx'^pj^^i^{x)dx. (7) 
Jo 

This is the mass that enters in the last term of equa- 
tion II]). Since the density profiles we consider do 
not have an outer edge where Phaio = 0, we truncate 
the cluster background halo at a maximum halo radius 
-^haio = 5 Mpc. Equation ([7]) is then solved numerically, 
to build an interpolation table for r in the range [0, R 
that is then used by the code. 

2.4. The Suhgrid Physics 

As mentioned in ^ 31. 31 there can be six possible phys- 
ical outcomes for our simulated cluster galaxies. In the 
following subsections, we describe the associated subgrid 
physics for each mechanism we use in our simulations. 
The possible outcomes are: (1) the galaxy merges with 
another galaxy ( ii2.4.ip . (2) the galaxy is destroyed by 
the tidal field of a larger galaxy but the fragments ac- 
crete onto that larger galaxy ( §2.4.4p . (3) the galaxy is 
destroyed by tides of a larger galaxy and the fragments 
are dispersed in the intracluster medium ( g2.4.3p . (4) the 
galaxy is destroyed by the tidal field of the background 
halo ( ti2.4.3p . (5) the galaxy survives to the present (i.e., 
it is not destroyed by any process), and (6) the galaxy 
is ejected from the cluster f §2.4.5p . We describe our ap- 
proach of simulating galaxy harassment in §2.4.21 

2.4.1. Encounter: Merger 

We simulate a pair of galaxies colliding (or synony- 
mously, having an encounter) and the further conse- 
quences (e.g., merging) in the following way. An en- 
counter is accounted for when two galaxies i and j, of 
radii Si and Sj, overlap such that the center of the galaxy 
j is inside the galaxy i. Numerically the criterion is 



< Si, where, 



Vj\ is the distance between 



the centers of the galaxies. If and Vj are the velocities 
of galaxies i and j, the center of mass velocity of the pair 
is Vein = {rriiVi + mj-Vj) / {nii + ruj). The kinetic energy 
in the center-of-mass rest frame is 



K, 



1 



1 I 



The gravitational potential energy of the pair is 



Gmimj 



(8) 



(9) 
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Even though we are treating each galaxy as a single 
particle, in reality a galaxy is a gravitationally bound 
system with an internal kinetic energy and a potential 
energy, and these energies must be included in the total 
energy of the interacting pair of galaxies. Considering a 
galaxy as a bound virialized system its internal energy is 



u, = a 



potential 



f/ki 



c/, 



kinetic 



potential 



CGmj 
2s, 



, (10) 



where C is a geometrical factor which depends on the 
mass distribution in the galaxies. Throughout this paper, 
we assume C = 1 (see Appendix A). 

We then compute the total energy of the galaxy pair 
(in the center of mass frame) as 

E,^^K,j+W,,+U, + U,. (11) 

If Eij < 0, i.e., the system is bound, we then allow 
the galaxies to merge to form a single galaxy of mass 
Wjnorgc = rrii + rrij. To compute its radius, we assume 
that energy is conserved, hence the total energy Eij in 
the center-of-mass rest frame is all converted into the in- 
ternal energy of the merger remnant. Its radius is then 
computed from equation (fTO|) . 



-^merged 



2 

merge 



CG(m, 



2U„ 



2\E,, 



(12) 

' merge ^ | J^i] | 

The position and velocity of this merged object are set 
to those of the center-of-mass values of the galaxy pair 
before merger, in order to conserve momentum. 

2.4.2. Encounter: Galaxy Harassment 

In a high-speed encounter the two interacting galaxies 
come into contact for a brief amount of time. The galax- 
ies might survive a merger or tidal disruption, but the 
encounter adds some internal energy into them, making 
them less bound. We refer to this process as Galaxy Ha- 
rassment. This process has been originally suggested as 
a possible explanation for t he origin of the mo rphology- 
density relation in clusters (jMoore "eTallflQQei) . We in- 
corporate galaxy harassment in our algorithm by increas- 
ing the radius (or the internal energy) of a galaxy when 
it experiences a non-merger encounter. This enlargement 
makes a galaxy more prone to tidal disruption at a next 
encounter. 

In equation pT|) . if Eij > 0, i.e., the system is not 
bound, then the galaxies in our simulation do not merge 
in the collision. Rather a part of the kinetic energy of 
the galaxies is converted into internal energy, making the 
collision inelastic. We assume that an equal amount of 
energy is transfered to each galaxy. Denoting the energy 
transfered as AE, the kinetic energy of the pair decreases 
by AE, while the internal energy of each galaxy increases 
by AE/2. We assume 



AE 

— =^ram{m,\U,\ 



(13) 



where rj is an energy transfer efficiency whose value is 
taken as ?7 = 0.2. The internal energies of the two galax- 
ies after the encounter are Uf^"' = [/before _^ ^^/2 and 
[/after ^ [/before ^_^^/2, respectively. By choosing < 1, 
we are ensuring that the internal energy of each galaxy 
remains negative, that is, the transfer of energy does not 
unbind the galaxies. We recalculate the velocities and 



Fig. 1. — Calculation of the effects of tides caused by a galaxy 
of mass nii on a galaxy of mass nij and radius sj . The two largest 
arrows show the gravitational accelerations caused by galaxy i; the 
two smallest arrows show the accelerations caused by galaxy j. See 
j2X3l for details. 



Vj after collision while conserving momentum, assuming 
that only the magnitudes of velocity change, the direc- 
tions remaining the same. We recalculate the size of each 
galaxy in the pair using equation (|12p . 



^ after 



CGm| 



^ after 



2\U' 



aftc 



2\U: 



afte 



(14) 



While allowing a size increase of the galaxies accord- 
ing to equations and p^ . we also considered a size 
cutoff. We assumed that the galaxies could grow only up 
to a maximum size given by the size of the largest galaxy 
at the beginning of the simulation. 

2.4.3. Tidal Disruption: Intracluster Stars 

We consider two possible sources of external gravita- 
tion for the tidal disruption of a galaxy j: another galaxy 
i, or the background cluster halo. The tidal force on a 
galaxy due to the gravitational field of the external source 
is meaningful only if the galaxy lies entirely on one side of 
the external source, when the tides are directed radially 
outwards tending to tear apart the galaxy. Our calcula- 
tion of the tidal field caused by a galaxy i of mass nii 
is illustrated in Figure [T] The galaxies i and j are sepa- 
rated by a distance rij . We calculate the resultant fields 
between two diametrically opposite points inside galaxy 
j, located at a radial distance d < Sj along the line join- 
ing the centers of the two galaxies. The two small and 
two large arrows in Figure [1] indicate the gravitational 
field (or acceleration) at the opposite points caused by 
galaxy j (self-gravity) and by galaxy i (external source 
of gravitation) , respectively. The magnitude of the tidal 
field is given by the difference between the gravitational 
field caused by galaxy i at the two opposite points,^ 



Grrii 



Grrii 



a: 



galaxy _ 

(ry - df (ry + df 



(15) 



The gravitational field caused by galaxy j (two small 
arrows in Figure [1]) is directed radially inwards and acts 
opposite to the tides, tending to keep the galactic mass 
inside radius d intact. The difference between that self- 
gravitational field at the two opposite points is 



2Gmj{d) 
d2 ■ 



(16) 



This reduces to the well-known form a: 



galaxy 
tide 



oc d/rf^ in the 



limit d <^ rij , but we do not make this approximation here. 
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where Tnj{d) is the mass of galaxy j inside radius d. 
When a^^^^^^ = Ograv then the tides will exceed self- 
gravity at radii larger than d, while self-gravity will ex- 
ceed the tides at smaller radii. Thus the layers of galac- 
tic mass located between radii d and sj would become 
unbound, while the ones located inside radius d would 
remain bound. Hence, the galaxy would be partly dis- 
rupted. In our code, we simplify things by using an 
"all-or-nothing" approach. A galaxy is either totally dis- 
rupted, or not disrupted at all. We consider that a galaxy 
is disrupted if half of its mass or more becomes unbound. 
If we assume an isothermal sphere density profile (as 
in Appendix A), then the half-mass radius is given by 
d = Sj/2. This is the value of d we use in equations (I15|) 
and (|16p . The criterion of tidal destruction then becomes 
<rde^^(d)>ag_(d),withd = 5,/2. 

We also consider tidal disruption by the background 
cluster halo, but only if rj > Sj (that is, the galaxy does 
not overlap with the center of the halo) . The magnitude 
of tidal field due to the cluster halo is 



halo 
hide 



GMhaio(rj - d) GA/haio(rj + d) 



(r, -f df 



(17) 



where Mhaio(?') is given by equation If ci}^^l°{d) > 
agrav('^), with d = Sj/2, galaxy j is tidally destroyed by 
the gravitational field of the halo. 

When galaxy j is considered to have been tidally de- 
stroyed by another galaxy i, the fragments of the dis- 
rupted galaxy might accrete onto galaxy i { W2AA\ . or 
they might be dispersed into the ICM when Eij > 
[Eij being the total energy of the galaxy pair given by 
eq. pT|) ]. For tidal destruction by the cluster halo the 
disrupted fragments are always dispersed into the ICM. 
In both cases, the destroyed galaxy is removed from the 
list of existing particles. The code keeps track of the 
amount of mass added to the ICM (in the form of IC 
stars) by tidal disruption. This quantity is initialized to 
zero at the beginning of the simulation, and every time 
a galaxy is tidally destroyed with its fragment dispersed, 
the mass of that galaxy is added up to the mass of IC 
stars. 

2.4.4. Tidal Disruption: Accretion 

We consider a possibility of accretion of the fragments 
of a tidally disrupted galaxy onto the galaxy causing the 
tides. This happens for the case of tidal disruption due 
to galaxies only (if the disruption is caused by the back- 
ground cluster halo, the fragments are always dispersed 
as IC stars). This situation occurs when the conditions 



a°; 



galaxy 



> a^^^^ and Eij < are both satisfied (see ^2.4.3p . 
The tidally disrupted galaxy accretes onto the more mas- 
sive galaxy. The mass of the bigger galaxy increases from 
TOi to rrii + nij . Thus a tidal disruption followed by ac- 
cretion is similar to a merger f ^2.4.ip . but these events 
are counted separately. 

2.4.5. Ejection 

When a galaxy ventures at distances larger than the 
cluster halo truncation radius R^^^ (see ^2.3^ . we con- 
sider that this galaxy has escaped from the cluster, and 
we remove it from the list. If we kept that galaxy, it 
might eventually return to the cluster. But in reality the 



universe contains many clusters, and a galaxy that moves 
sufficiently far away from one cluster will eventually feel 
the gravitational influence of other clusters, something 
that our algorithm, which simulates an isolated cluster, 
does not take into account. As we shall see, the ejec- 
tion of galaxies from a cluster is quite uncommon in our 
simulations. 

3. THE SIMULATIONS 
3.1. Cosmological Model 

We consider a ACDM model with the present matter 
density parameter, ^Im — 0.241, baryon density param- 
eter, ilb = 0.0416, cosmological constant, J^a = 0.759, 
Hubble constant, = 73.2kms~^Mpc"^ [h = 0.732), 
primordial tilt, Us = 0.958, and CMB temperature, 
TcMB = 2.725 K, consistent with the results of WMAP3 
(Spe rgel et al.l [20071. Even though the simulations pre- 
sented in this paper are not "cosmological" (we simulate 
isolated, virialized clusters), the particular choice of cos- 
mological model enters the picture twice: in the deter- 
mination of the radii of galaxies (see next section), and 
in the calculation of the elapsed time between the initial 
and final redshifts of the simulation. 

In each simulation a cluster is evolved from z = 1 to 
the present [z — 0). We assume that the cluster will 
not experience any major merger during this period, and 
that, therefore, it is a good approximation to treat it as 
isolated. For our ACDM model, this represents a total 
evolutionary time of 7.63 Gyr. 

3.2. Initial Conditions 

To set the initial conditions of our simulations, we need 
to determine the initial mass m, radius s, position r, 
and velocity v of each galaxy. To determine the mass, 
we first assume that the luminosities of galaxies are dis- 
tributed accordi ng to the Schechter luminosity function 
(ISchechtedH97l . 



4>{L)dL = r — 



-L/L'^ 



(18) 



Here we use L* = 3.097 x 1G^°Lq (corresponding to 
absolute magnitude M^^ = -20.07), and a = -1.28, 
which is appro priate for galaxies located in clusters 
()De Propris et a l. 2003) . This luminosity function spans 
over —22.5 < Mb, < —15. While it might be rea- 
sonable to assume fixed values of L* and a, the value 
of (j)* most likely varies amongst clusters. So we nor- 
malize equation (|18p by imposing that, in each cluster, 
there are A^o galaxies with luminosities L > Lq. We use 
No = 25. and Lp = 0.2L* (corresponding to Mb = -19) 
(| Lewis et al.l |2002|) . We select the luminosities using 
a Monte Carlo rejection method. We then assume a 
constant m ass-to-light ratio T = 193 h Mq/Lq = 73 
(jBrainerd fc Specian 20Q3|), and convert the luminosities 
to masses. The Schechter function spans up to a max- 
imum mass Mmax = 220 x IO^^Mq. To generate the 
dwarf galaxies, the same Schechter function is extrapo- 
lated up to a minimum mass Mj^m = 1 x 1O^M0. 

We take the radius s of each galaxy to be equal to 
the virial radius r2oo (radius containing matter with 200 
times the mean density of the Universe at the epoch of 
galaxy formation) corresponding to the galaxy mass m = 



Fig. 2. — Initial conditions for run A12. Top panel: initial conditions at ^ = 1. The solid circles indicate the virial radii of galaxies. 
The large circle is the maximum distance r = 3Ro = 2.08 Mpc from the cluster center. Lower left panel: same as top panel, with symbols 
rescaled to optical diameter of real galaxies. Bottom right panel: enlargement of the central (0.6 Mpc)^, (box on lower-left panel). 



M200 using 



M: 



SOOtt 



200 



''200P(1 + ^coll) 



(19) 



Here, p = Pait^M = SH^SIm /^ttG is the mean matter 
density in the present universe, and ZcoU is the redshift of 
collapse when the galaxy formed. To obtain Zcoib we use 
a simple spherical collapse model. First, by filtering the 
power spectrum for our ACDM model, we calculate the 
standard deviation a{m) of the linear density contrast 
S = [p — p)/p at the mass scale m. The distribution of 
the values of 5 is then given by a Gaussian, 



V[5) cx exp ( - — 



(20) 



We pick randomly a present density contrast 5q — 6{z = 
0) from this distribution, using a Monte Carlo rejection 
method, and solve the following equation to get the col- 
lapse redshift ZcoU, 



A. 



^+(^:cq1i) 

' '5+(0) 



(21) 



where (^+(2;) is the linear growing mode (for A 7^ mod- 
els, see, e.g.. lMartelllT99l . Here, = 3(127r)2/3/20 = 
1.686 is the overdensity predicted by linear theory at rec- 
oUapse. We solve this equation numerically for Zcoib ^-ud 



substitute the solution in equation (jl9p . which we then 
solve to get the radius s = r2oo ■ 

To determine the locations of galaxies inside a cluster, 
we assume that their distribution is isotropic (in a statis- 
tical sense). We can therefore choose the spherical coor- 
dinates {9, 4>) of each galaxy randomly, using ip = 27rX0, 
and cos 9 — 2Xg — 1 , where and Xg are random num- 
bers drawn from a uniform deviate between and 1. We 
still need to determine t he radial coordinate r . Using 
the CNOC cluster survey, iCarlberg et all (|1997D showed 
that the radial mass density p(r) of matter and the radial 
number density iy{r) of galaxies are roughly proportional 
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to each other, where b oth p(r) and t ^(r) a re approxi- 
mated by NFW profiles. iGirardi et all found that 
the halo mass follows the galaxy distribution in clusters, 
using a /?-model for the halo/galaxy volume density pro- 
file. We assume that this proportionality holds for all 
clusters, and we generalize it to to all the density profiles 
we use. Thus, the assumed background halo mass den- 
sity profile [eqs. ([3]) and (|4])] gives us v{r). We can then 
select the initial distances r bom the cluster center using 
again a Monte Carlo rejection method. Since the masses 
and locations of the galaxies have been determined sep- 
arately, we need to pair them, that is, for each selected 
location, to decide which galaxy goes there. We do ex- 
pect the most massive galaxies to reside near the center 
of the cluster. However, low-mass galaxies are not all lo- 
cated at large radii, and some of them might be located 
in the central region of the cluster as well. Indeed, if the 
galaxies in the central region were all massive, it would 
be impossible to reproduce the desired number density 
profile i'{r) and not have the galaxies overlap. 

To prevent any overlap, we locate the galaxies as fol- 
lows: we first position the most massive galaxy at the 
center of the cluster. Then we locate the next 7 most 
massive galaxies between radii Rq and 2Ro, where Rq 
is 3 times the radius of the most massive galaxy. We 
then locate the next 19 most massive galaxies between 
radii 2i?o and 3i?o- Finally, the remaining galaxies are 
located randomly between radii and 3i?o- During the 
process we check that the galaxies do not overlap, by en- 
suring that the distance between the edges of two galax- 
ies is greater than the radius of the larger galaxy, i.e., 
fij — Si — Sj > max(si, Sj). In the process of locating a 
galaxy, if this criterion is not satisfied, we simply reject 
that location and generate a new one. 

After assigning the masses, radii, and positions of all 
the galaxies, we determine the velocity of each galaxy. 
We consider the velocity a galaxy would have if it were 
in a perfect circular orbit at radius r, 

1/2 

G ' 



- I A4alo(»^) 

r ' 



(22) 



where the sum only includes galaxies inside radius r. The 
norm of the velocity is chosen by giving a random 10% de- 
viation to the circular velocity, i.e., v = fcirc(l + 0.1X„), 
where Xy is a random number between —1 and 1. For 
the direction of the velocity, we follow a similar random 
angle generation technique as we did for the positions of 
galaxies. 

Figure [2] illustrates the initial conditions for one of our 
simulations (run A12). The top panel shows the cluster 
at z = 1. The large circle represents the maximum dis- 
tance 3i?o within which the galaxies are located initially. 
Each dot represents a galaxy, with the most massive one 
located in the center. Even though massive galaxies tend 
to be larger, there is no direct correspondence between 
the masses and radii because of the dependence on Zcou 
in equation (|19p . whose determination involves a Monte 
Carlo method. 

Visually, this looks quite different from the optical pic- 
tures of actual clusters like Virgo. This is because each 
dot has a radius s equal to the virial radius r200 7 that 
can exceed the optical radius by an order of magnitude. 
In the bottom left panel of Figure [2l we show the same 



cluster, with all the dots rescaled in size so that the an- 
gular diameter of the central galaxy is equal to 8.3' at a 
distance of 16.8 Mpc, which is the observed optical diam- 
eter of M87. The bottom right panel shows a zoom-in of 
the central cluster region. It looks qualitatively similar 
to pictures of the central region of Virgo. 

4. RESULTS 

We started by performing 9 series of simulations, for a 
total of 132 simulations. Table 1 summarizes the char- 
acteristics of each series. The first 2 columns show the 
series name and the number of runs, respectively. The 
slope of the Schechter luminosity function at z = 1 (used 
to generate the initial conditions) is listed in column 
3. Columns 4 to 8 give the characteristics and rele- 
vant parameter values of the background cluster halo 
profile. Columns 9 and 10 indicate respectively whether 
a cD galaxy was included in the cluster simulation, and 
whether galaxy harassment was included as part of the 
subgrid physics. 

4.1. Series A: Initial Simulations 

We performed an initial series of 16 simulations, us- 
ing for the background halo a /9-profile with (3 — 0.33, 
a core radius Tc = 3kpc, and a central density po = 
8.14 X 10~ ^^gcm~'^, which is appropriate for a cluster 
like Virgo (|Piffaretti fc Kaas tra 2006). For this series, 
we did not include galaxy harassment. Our results are 
shown in Table [21 It shows the run number (column 1), 
and, at the beginning of the run, the total mass Mtotai in 
galaxies, in units of lO^^M© (column 2), the number of 
galaxies A^totai (column 3), and the Schechter luminosity 
function [eq. (jlSp ] exponent ctsta.rt (column 11). This ex- 
ponent was obtained by performing a numerical fit to the 
distribution of galaxy masses. Because the masses were 
determined from a Monte Carlo rejection method, the 
exponent can differ from the intended value a = —1.28 
in equation (fTS]) and listed in Table [l] but the deviations 
are small. Averaging over all runs, we get 

= -1.280 ± 0.020. (23) 

Columns 4 — 8 in Table [2] show the number of galax- 
ies N„ 



N, 



gal 
tides 



destroyed by mergers, the number of galaxies 

destroyed by tides caused by a massive galaxy, 
with the fragments dispersed in the ICM, the number 
of galaxies A^accr destroyed by tides caused by a mas- 
sive galaxy, with the fragments being accreted onto that 
galaxy, the number of galaxies A^t^dos destroyed by the 
tidal field of the background halo, and the number of 
galaxies A'ojcct ejected from the cluster, respectively. Col- 
umn 9 shows the fraction by numbers of galaxies /surv 
that survive to the present. 

We did not find a single occurrence of a galaxy de- 
stroyed by tides from the background halo, and the num- 
ber of galaxies ejected is either or 1. There are large 
variations in the other numbers from one run to the next, 
but some trends are apparent. Typically, 50% to 60% of 
the galaxies are destroyed. Run A2 is an extreme case, 
with 78% of the galaxies being destroyed. The destruc- 
tion of galaxies by mergers dominates over the destruc- 
tion by tides, by more than a factor of 2 except for run 
A7. If we treat the cases of tidal disruption followed by 
accretion as being mergers, then mergers dominate even 
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TABLE 1 

Series of Simulations 



Series 


Runs 


'^start 


profile 


13 


Po.Ps [gem ^] 


c 


rc , [kpc] 


cD 


Harassment 


A 


16 


-1.28 


/3- Virgo 


0.33 


8.14 X 10"^^ 




3 


X 


X 


B 


17 


-1.28 


/3- Virgo 


0.33 


8.14 X 10-26 




3 


X 


V 


C 


17 


-1.36 


/3- Virgo 


0.33 


8.14 X 10-26 




3 


X 


V 


D 


16 


-1.36 


/3- Virgo 


0.33 


8.14 X 10-26 




3 


V 


V 


E 


16 


-1.36 


/3-Perseus 


0.53 


7.27 X 10-^6 




28 


X 


V 


F 


16 


-1.36 


/3-Perseus 


0.53 


7.27 X 10-26 




28 




V 


G 


10 


-1.28 


NFW 




2.35 X 10-2^' 


5 


200 


X 


V 


H 


14 


-1.31 


NFW 




2.35 X 10-25 


5 


200 


X 


V 


I 


10 


-1.31 


NFW 




2.35 X 10-25 


5 


200 


V 


V 



TABLE 2 

Simulations for Series A 



Run 


Mtotal [10" Mq] 


-'Vtotal 


^ merge 


tides 


^ accr 


^halo 
tides 


A'cjcct 


f surv 


Acs 


'-'start 


end 


Al 


855.0 


440 


137 


54 


6 





1 


0.550 


0.225 


-1.26 


-1.23 


A2 


862.0 


702 


343 


182 


24 





1 


0.217 


0.490 


-1.28 


-1.12 


A3 


1100.4 


480 


126 


43 


7 








0.633 


0.136 


-1.28 


-1.24 


A4 


831.0 


530 


165 


67 


11 








0.542 


0.148 


-1.26 


-1.15 


A5 


980.4 


459 


175 


58 


8 





1 


0.473 


0.128 


-1.25 


-1.18 


A6 


967.8 


640 


230 


74 


5 








0.517 


0.157 


-1.29 


-1.23 


A7 


720.0 


579 


263 


114 


7 








0.337 


0.359 


-1.28 


-1.20 


A8 


757.8 


435 


200 


66 


51 





1 


0.269 


0.320 


-1.27 


-1.26 


A9 


925.6 


514 


223 


83 


12 





1 


0.379 


0.245 


-1.28 


-1.17 


AlO 


858.9 


525 


204 


54 


8 








0.493 


0.257 


-1.31 


-1.24 


All 


880.2 


547 


224 


68 


9 





1 


0.448 


0.257 


-1.33 


-1.25 


A12 


826.1 


548 


255 


72 


15 





1 


0.374 


0.225 


-1.29 


-1.18 


A13 


1041.5 


431 


166 


51 


6 








0.483 


0.275 


-1.27 


-1.25 


A14 


957.3 


486 


174 


52 


4 








0.527 


0.261 


-1.29 


-1.22 


A15 


860.5 


520 


160 


64 


15 





1 


0.538 


0.260 


-1.26 


-1.18 


A16 


858.0 


483 


255 


72 


8 








0.306 


0.323 


-1.28 


-1.19 



more over tidal disruption. When galaxies are destroyed 
by tides, the dispersion of fragments into the ICM al- 
ways dominates over the accretion of fragments onto the 
massive galaxy, but the ratio varies wildly, from 114:7 in 
run A7 to 66:51 in run A8. 

We can now evaluate the fraction /j^g of the total lu- 
minosity of the cluster that comes from intracluster stars. 
Since we assume a constant mass-to-light ratio, this frac- 
tion is given by 



/ics — 



71 ,rgal I /urhalo 



M. 



total 



M, 



(24) 



eject 



where the letter M refers to the mass in galaxies, rather 
than their number. The galactic mass contribution to 
the ICM consists of galaxies destroyed by tides of another 
more massive galaxy, and by tides of the background halo 
(though there are no such cases in this series). 

The values of /j^g are listed in column 10 of Table O 
Again, there are large variations. In particular, the frac- 
tion is very large for run A2, and very small for runs A3 
and A5. Averaging over all runs, we get 



/ics 0.254 ±0.093. 



(25) 



Even though, in most cases about half the number of 
galaxies are destroyed, they tend to be low-mass galaxies, 
which explains why /ics < 1 — /surv, for all the runs. 

The galaxies being destroyed by mergers and tides, or 
escaping are mostly low-mass galaxies. This leads to a 
flattening of the galaxy mass distribution function. We 



computed the best numerical fit to the Schechter lumi- 
nosity function exponent a [eq. (|18p ] for the surviving 
galaxies at the end of the simulations. This is listed as 
a^jjjj in column 12 of Table [2l Averaging over all runs. 



we get 



a^^^j = -1.206 ±0.040. 
4.2. Series B: Turning on Harassment 



(26) 



We modified the algorithm to include the effect of 
galaxy harassment (see !j2.4.2p . and rerun the calcula- 
tions of series A with the same initial conditions. We 
also added one more run, B17. The results are shown in 
Table 131 which follows the same format as Tabled Com- 
paring with series A, the number of galaxies destroyed by 
mergers is very similar, but the number of galaxies de- 
stroyed by tides tends to be significantly higher. For in- 
stance, it goes from 67 to 94 for runs A4-B4, and from 64 
to 86 for runs A15-B15. This is because, when a galaxy is 
subjected to harassment, its binding energy is reduced, 
and it becomes more prone to experience tidal disrup- 
tion later. The number of tidal disruptions followed by 
accretion does not change significantly, though. Hence, 
the additional, tidally-disrupted galaxies almost all con- 
tribute to the intracluster stars. The values of /ics are 
therefore increased relative to series A. The mean value 



IS 



/] 



ICS 



0.269 ±0.081 



(27) 



This is not significantly larger than for series A. The ad- 
ditional galaxies destroyed are mostly low-mass galaxies. 
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TABLE 3 

Simulations for Series B 



Run 


A^total[lOl^M0] 


-'Vtotal 


^ merge 


I 

tides 


^ aecr 


tides 


-^eject 


f surv 


Acs 


'-'start 


end 


Bl 


855.0 


440 


134 


62 


g 





I 


0.539 


0.237 


— 1.26 


— 1.24 


B2 


862.0 


702 


337 


210 


25 





I 


0.184 


0.454 


— 1.28 


— 1.11 


B3 


1100.4 


480 


126 


48 


y 


Q 


Q 


0.623 


0.139 


— 1.28 


— 1.24 


B4 


831.0 


530 


159 


94 


12 


Q 


Q 


0.500 


0.184 


— 1.26 


— 1.15 


B5 


980.4 


459 


171 


78 


7 





I 


0.440 


0.154 


— 1.25 


— 1.21 


B6 


967.8 


640 


226 


85 


g 


Q 


Q 


0.505 


0.173 


— 1.29 


— 1.23 


B7 


720.0 


579 


249 


134 


9 








0.323 


0.364 


-1.28 


-1.20 


B8 


757.8 


435 


172 


79 


62 





1 


0.278 


0.271 


-1.27 


-1.27 


B9 


925.6 


514 


223 


100 


12 





1 


0.346 


0.270 


-1.28 


-1.19 


BIO 


858.9 


525 


206 


68 


7 








0.465 


0.298 


-1.31 


-1.25 


Bll 


880.2 


547 


233 


85 


9 





1 


0.400 


0.280 


-1.33 


-1.27 


B12 


826.1 


548 


241 


92 


19 





1 


0.356 


0.240 


-1.29 


-1.20 


BIS 


1041.5 


431 


160 


58 


4 








0.485 


0.269 


-1.27 


-1.26 


B14 


957.3 


486 


173 


59 


5 








0.512 


0.274 


-1.29 


-1.22 


B15 


860.5 


520 


155 


86 


17 





1 


0.502 


0.277 


-1.26 


-1.19 


B16 


858.0 


483 


264 


90 


6 








0.255 


0.354 


-1.28 


-1.16 


B17 


792.6 


451 


178 


83 


16 





1 


0.384 


0.346 


-1.24 


-1.13 



We also recalculated the best-fit Schechter exponent a 
for the surviving galaxies at z — 0. The mean value for 
the runs in this series is 

a^nd = -1-207 ±0.048. (28) 

Figure [3] shows the total galaxy counts in mass bins, 
obtained by combining all the runs in series B, along 
with the fitting curves to a Schechter distribution func- 
tion [eq. (fT8| ]. The best fit Schechter exponent for the 
initial galaxy distribution (the upper curve at z = 1) is 
a = —1.28, and for the final surviving galaxy distribution 
(the lower curve at z = 0) is a = —1.20. These values of 
a were obtained by performing the numerical Schechter 
function fits on the combined set of galaxies taken from 
all the 17 runs in this series, which amounts to 8770 
initial galaxies at z = 1 and 3614 surviving galaxies at 
z = 0. 

Clearly from Figure [3l the fit at z = (lower curve) 
is excellent. This shows that, in our simulations, a 
Schechter mass (luminosity) distribution function at z = 
1 remains a Schechter distribution all the way to z = 0, 
though half of the galaxies are destroyed. Only the slope 
a changes with time. 

4.3. Series C: Steeping up the Mass Distribution 
Function 

In the simulations of series A and B, the Schechter ex- 
ponent a evolves from a ~ —1.28 at z = 1 to a ~ —1.21 
at z = 0. Analyzing the combined set of galaxies of all 
the 17 runs in series B, we obtained the best fit Schechter 
exponent for the surviving galaxy distribution at z = 
as a = —1.20. This is a problem, since the value of 
a = -1.28 is based on . observations of nearby clusters 
(jPe Propris et all 120031 ) , and should be valid for clus- 
ters of galaxies at z = 0, whereas we used the same 
a = —1.28 to start our simul ations at z = 1, an d it fiat- 
tened to a = -1.20 at z = 0. iRvan et aP (|2007t l recently 
determined the luminosity function of a large sample of 
galaxies at z ~ 1, and concluded that there is a steepen- 
ing of the faint-end slope with redshift, which is expected 
in the hierarchical formation scenario of galaxies. They 
obtained a value of the faint-end slope a = —1.32 ± 0.07 
at z = 1. 



1400 




Log,o [Galaxy Mass (10" U^„„)] 

Fig. 3. — Mass distribution function for galaxies in Series B, 
obtained by adding the numbers for all runs. Results are plotted 
for initial 8770 galaxies at 2: = 1 (asterisks) and surviving 3614 
galaxies at 2 = (plus signs). The cu rves show the best fit of a 
Schechter distribution function (eg. I18| l. with « = — 1.28 at 2 = 1 
(upper curve), and o = — 1.20 at 2 = (lower curve). 

In our simulations we take into account this fiattening 
of the luminosity function over time as explained below. 
The results of series A and B suggest that |a| decreases 
by ^ 0.08 between z = 1 and z = 0. Using this as a 
guide, we performed a new series of simulations, series C, 
using ttstart = -1-36, with the hope that this value will 
evolve toward something close to a ~ —1.28 at z = 0. 
The results are shown in Table H) The average values of 
a are 

a,tart = -1-357 ±0.021, (29) 
a„,d = -1-272 ±0.050. (30) 

Figure [4] (analogous to Figure [3]) shows that at z = 0, 
a Schechter distribution function is still a good approxi- 
mation to the mass distribution. Combining the galaxies 
from all the 17 runs in this series (13628 initial galaxies 
at z = 1, and 5356 surviving galaxies at z = 0), the best 
fit Schechter exponent for the initial galaxy distribution 
(the upper curve) is a = —1.36, and for the surviving 
galaxy distribution (the lower curve) is a = —1.27. This 
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TABLE 4 

Simulations for Series C 



Run 


A^total[lO"M0] 


.'Vtotal 


^ merge 


I 

tides 


^ aecr 


tides 


.^eject 


f surv 


Acs 


'-'start 


end 


CI 


1144.3 


948 


322 


123 


13 





7 


0.509 


0.162 


— 1.34 


— 1.27 


C2 


1009.9 


936 


283 


138 


19 





1 


0.529 


0.195 


— 1.38 


— 1.30 


C3 


914.2 


603 


316 


97 


13 


Q 


Q 


0.294 


0.265 


— 1.35 


— 1.24 


C4 


970.5 


635 


162 


86 


11 


Q 


1 


0.591 


0.216 


— 1.32 


— 1.28 


C5 


935.9 


672 


250 


101 


7 





1 


0.466 


0.158 


— 1.33 


— 1.25 


C6 


948.2 


1136 


502 


268 


g 


Q 


Q 


0.317 


0.318 


— 1.35 


— 1.21 


C7 


926.4 


884 


276 


155 


10 





1 


0.500 


0.258 


-1.35 


-1.28 


C8 


818.4 


696 


347 


142 


62 





1 


0.207 


0.423 


-1.35 


-1.26 


C9 


875.9 


903 


421 


255 


44 





1 


0.202 


0.405 


-1.34 


-1.19 


CIO 


1110.5 


785 


286 


100 


20 





1 


0.482 


0.168 


-1.39 


-1.39 


Cll 


864.2 


809 


338 


139 


18 





1 


0.387 


0.386 


-1.40 


-1.36 


C12 


1053.2 


978 


395 


212 


29 








0.350 


0.294 


-1.36 


-1.29 


C13 


1151.7 


790 


344 


127 


24 





1 


0.372 


0.262 


-1.37 


-1.29 


C14 


876.8 


684 


335 


133 


23 





1 


0.281 


0.320 


-1.37 


-1.21 


C15 


788.0 


769 


240 


140 


27 








0.471 


0.243 


-1.35 


-1.26 


C16 


905.1 


578 


284 


80 


24 





1 


0.327 


0.430 


-1.36 


-1.29 


C17 


889.2 


822 


325 


169 


13 








0.383 


0.323 


-1.36 


-1.25 



value of a is close enough to our target value of —1.28. 
So from now on, in all subsequent series with the (3 model 
halo density profile, we will use an initial a of —1.36, as 
shown in T able [U Thi s valu e is well inside the range 
obtained bv lRvan et"aI1 ()2007D . 

Using a steeper galaxy distribution while still requiring 
that the clusters contain 25 galaxies with L > 0.2L* (see 
initial conditions in results in the initial number of 
galaxies being larger by a factor of about 2 (column 3 
of Table [4|) . But the numbers of galaxies destroyed by 
mergers and tides are also higher relative to series B. As 
a result the trends are similar. In particular, mergers 
still dominate over tides by more than a factor of 2. 

The run CI has a larger number of galaxies ejected 
from the cluster. This is because the most massive 
galaxy, located at the center of the cluster, was partic- 
ularly large. Its radius was s = 385 kpc, compared to 
s < 300 kpc for the other runs. This increased the value 
of Rq (see used for setting up the initial conditions. 
As a result, more galaxies were located at larger radii, 
where they are more likely to escape. 

The mean value of /jcs for this series is 

/ics = 0.284 ±0.090. (31) 

4.4. Series D: Adding a cD Galaxy 

A cD (central dominant) galaxy is a very bright super- 
giant elliptical galaxy with an extended env elope (or a 
diffus e halo) found at the center of a cluster (jSchombertj 
I1988D . Several galaxy clusters have been found to have 
cD g al a.xies at their c ei iters (e. g., [Q uintana & Lawrie 
198a iHill fc Oegerld [l998l: lOegerle fc Hil| l200ll : 
Jordan et al.l l2004t ISeigar. Graham, fc JerienI I2007D . 
We performed some simulations by incorporating a 
cD galaxy in the clusters. Being the brightest and 
most massive cluster galaxy, the mass of a cD is 
larger than the prediction of the normal Schechter 
distribution function (eq. I18p . So we introduced the 
cD galaxy manually at our simulated cluster center. 
We adopted a luminosity of LcD — lOL*, which is 
a canonical value for a cD. Using our mass to light 
ratio (T — 193 h Mq/Lq, ^3.2^ . this corresponds to 
a cD galaxy mass of M^b = 437.6 x 10" Mq. When 
we wanted a cD galaxy present in the simulation we 
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Log,„ [Galaxy Mass (10" M,J] 

Fig. 4. — Same as Figure [S] for series C. Results are plotted 
for initial 13628 galaxies at z = 1 (asterisks) and surviving 5356 
galaxies at z = (plus signs). The best fit Schechter distribution 
functions are with a = —1.36 at z = 1 (upper curve), and a = 
— 1.27 at 2 = (lower curve). 

changed the mass of the cluster central galaxy (see ij3.2p 
to the cD mass, M^u- This allowed us to keep the 
appropriate initial galaxy distribution for a cluster while 
incorporating a cD galaxy at rest, located at the center. 

We performed simulations by adding a cD galaxy to 
our Virgo-like cluster, and called it series D. The results 
are listed in Table [5l from which certain trends are clear 
after incorporating a cD galaxy in the simulation. The 
total galaxy mass increases since a massive cD galaxy 
is being added. More prominent than in the previous 
series A, B, and C, here galaxy mergers outnumber tides 
by factors ~ 2 — 3, which go as high as 4 in run D12. 

A striking new feature in cases incorporating a cD 
galaxy is the increase in the number of accretions af- 
ter tidal disruption by a galaxy, fully 1/4 of the galax- 
ies being acreted in run D5. Since in these accretions, 
the smaller galaxy is tidally destroyed and is absorbed 
(or merged) by the massive galaxy fi j2.4.4p . it appears, 
in our simulated clusters, that in the presence of a cD 
galaxy the number of effective mergers is very high. 

The mass fraction imparted to ICS has decreased in 
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TABLE 5 

Simulations for Series D 



Run 


A^total[lOl^M0] 


-'Vtotal 


^ merge 


I 

tides 


^ aecr 


tides 


-^eject 


f surv 


Acs 


'-'start 


end 


Dl 


1394.6 


1023 


345 


165 


83 








0.420 


0.127 


— 1.38 


— 1.33 


D2 


1374.4 


866 


330 


97 


37 





I 


0.463 


0.085 


— 1.38 


— 1.33 


D3 


1431.1 


906 


257 


199 


187 


Q 


Q 


0.290 


0.125 


— 1.39 


— 1.32 


D4 


1163.2 


754 


211 


152 


120 


Q 


Q 


0.359 


0.117 


— 1.33 


— 1.26 


D5 


1239.5 


788 


283 


167 


197 








0.179 


0.167 


— 1.33 


— 1.21 


D6 


1158.8 


527 


131 


59 


53 





1 


0.537 


0.065 


-1.34 


-1.30 


D7 


1300.7 


663 


236 


85 


73 





1 


0.404 


0.151 


-1.37 


-1.31 


D8 


1298.4 


773 


296 


113 


116 








0.321 


0.120 


-1.36 


-1.33 


D9 


1299.5 


652 


256 


95 


63 





1 


0.364 


0.141 


-1.39 


-1.33 


DIO 


1169.4 


589 


193 


76 


67 





1 


0.428 


0.106 


-1.37 


-1.32 


Dll 


1145.8 


741 


282 


143 


79 





1 


0.318 


0.124 


-1.38 


-1.25 


D12 


1227.5 


732 


235 


57 


22 





9 


0.559 


0.075 


-1.39 


-1.35 


D13 


1334.2 


952 


330 


166 


70 








0.406 


0.151 


-1.35 


-1.31 


D14 


1286.0 


844 


283 


184 


115 





1 


0.309 


0.203 


-1.34 


-1.24 


D15 


1470.8 


888 


298 


115 


56 








0.472 


0.140 


-1.36 


-1.33 


D16 


1454.4 


966 


288 


170 


125 








0.396 


0.086 


-1.34 


-1.30 



all the runs, with a value as low as 0.065 in run D6. 
To explain such a result, we note that the most massive 
central galaxy (cD or otherwise) in our simulated cluster 
is never destroyed because of its large mass. In an en- 
counter, it is normally the lower-mass galaxy that gets 
destroyed. Also the initial conditions of the most mas- 
sive galaxy (at rest at the center, see ^'S.2\i make it less 
likely to be destroyed by the tidal field of the halo. If the 
central galaxy is a cD, a large mass fraction (as high as 
38% in run Dll) is locked into it, which can never con- 
tribute to the ICS. So a smaller mass fraction is available 
to be transferred to the ICS, which eventually leads to a 
decrease in /ics- 
The mean value of /ics for series D is 

/ics = 0.124 ±0.036. (32) 

4.5. Series E & F: Other (3 Profiles 

In the next two series of runs, we consider a different 
background halo. We use a /3-profile with (i = 0.53, 
a core radius rc = 28kpc and a central density po = 
7.27 X 10~^^ gcm~'^, which is appropria te for a cluster 
like Perseus (jPiffaretti fc Kaastral [2006D . Series E and 
F do not include, and include a cD galaxy, respectively 
(hence series E should be compared with series C, and 
series F with series D). 

The results for series E are shown in Table [6] The 
most notable feature is that some galaxies are destroyed 
by the tidal field of the background halo, something that 
never happened with Virgo- like clusters. In order to ex- 
plain such a behavior we note that tidal disruption by 
the cluster halo generally occurs with galaxies at a dis- 
tance r <1 Mpc from the cluster center. Our simulated 
Perseus-like cluster halo mass profile rises more steeply 
than the Virgo-like cluster up to ~ 1.7 Mpc, making 
Perseus more massive in the inner regions. So a galaxy 
at a smaller distance, precisely at r < 1.7 Mpc, from 
the cluster center feels a larger tidal field from a more 
massive halo, and is more prone to be disrupted in the 
Perseus-like cluster. 

The numbers of other galaxy outcomes are compara- 
ble for Perseus-like and Virgo-like clusters, with mergers 
dominating over tides. The mean /ics for series E is 

/j(.g = 0.344 ± 0.093 . (33) 



This /ics is somewhat larger than the Virgo-like cluster 
mean (series C). This can be attributed to the non-zero 
tidal disruption by the cluster halo, resulting here in a 
finite contribution to the ICS mass fraction. 

Table [7] shows the results for series F, i.e., simulations 
of a Perseus-like cluster with a cD galaxy at the center. 
Here few galaxies are destroyed by the tidal field of the 
cluster halo, yet the number is smaller than in series E. 
It appears, then that the presence of a cD galaxy reduces 
the number of tidal disruptions by the background halo, 
since galaxies that would be destroyed by the tidal field 
of the central parts of the halo are being destroyed by 
the cD galaxy instead. 

Comparing the results for series D and series F (Virgo- 
like and Perseus-like clusters with a cD galaxy,) the 
numbers- merger, galaxy-tide and accretion are similar. 
Series F continues the trend of increased accretions when 
a cD galaxy is introduced. Also series F has a smaller 
fraction of mass going to ICS. The mean /ics for series F 
is 

/ics = 0.128 ±0.031 . (34) 

This /ics is very similar to that of the relevant Virgo-like 
cluster mean (series D). This implies that in the presence 
of a cD galaxy, the ICS mass fraction is not so sensitive 
to the parameters of the /3-model density profile. 

4.6. Series G, H, & I: NFW Profile 

We now consider a background halo described by a 
NFW profile [see N2.31 eq. ([U], with a scale radius 
rg = 200 kpc, and a concentration parameter c = 5. 
These values are adopted from obser vational studies 
of galaxy cluste rs (Arabadiis. B autz. fc Garmird 120021 : 
iPratt fc Arnaudl[2005. : ,Maughan et al.ll2007f) where the 
authors found the best fitting NFW model parameters 
for cluster mass profiles. 

We do not necessarily expect the flattening of the 
Schechter mass function to be the same for the NFW 
profile halo and the /3-profile halo. So at first we per- 
formed a series with a = —1.28 (see Table [1]), and called 
it series G. The results are listed in Table HI 

To contrast a NFW-model cluster with a /3-model clus- 
ter. Series G should be compared with Series B, since 
these are with a — —1.28, include galaxy harassment, 
and no cD galaxy. The most striking feature is the large 
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TABLE 6 

Simulations for Series E 



Run 


A^total[lO"M0] 


.'Vtotal 


^ merge 


I 

tides 


^ accr 


tides 


.^eject 


f surv 


Acs 


'-'start 


end 


El 


1065.4 


1023 


406 


159 


24 




I 


0.415 


0.297 


— 1 


.38 


— 1.31 


E2 


1034.0 


866 


326 


104 


7 


g 


2 


0.486 


0.168 


— 1 


.38 


— 1.31 


E3 


1051.3 


906 


405 


213 


19 


14 


Q 


0.282 


0.446 


— 1 


.38 


— 1.27 


E4 


799 4 


754 


250 


177 


35 




Q 


0.378 


0.391 


-1 


.34 


— 1.25 


E5 


868.7 


788 


394 


201 


42 


12 





0.176 


0.492 


-1 


.32 


— 1.22 


E6 


795.8 


527 


154 


67 


7 


6 


1 


0.554 


0.303 


-1 


.34 


-1.29 


E7 


947.7 


663 


276 


97 


9 


3 


1 


0.418 


0.376 


-1 


.37 


-1.31 


E8 


941.0 


773 


395 


121 


11 


8 





0.307 


0.348 


-1 


.37 


-1.28 


E9 


934.5 


652 


263 


109 


10 


6 


1 


0.403 


0.325 


-1 


.39 


-1.34 


ElO 


807.5 


589 


210 


98 


11 


5 


1 


0.325 


0.448 


-1 


.38 


-1.35 


Ell 


760.7 


741 


349 


142 


9 


8 


1 


0.313 


0.283 


-1 


.48 


-1.23 


E12 


946.4 


732 


227 


57 


4 


2 


28 


0.566 


0.149 


-1 


.39 


-1.35 


E13 


953.4 


952 


344 


180 


14 


11 





0.423 


0.349 


-1 


.35 


-1.29 


E14 


950.5 


844 


325 


197 


23 


6 


1 


0.346 


0.395 


-1 


.33 


-1.29 


E15 


1145.5 


888 


320 


118 


13 


5 





0.486 


0.339 


-1 


.36 


-1.29 


E16 


1125.4 


966 


439 


146 


10 


10 


1 


0.373 


0.393 


-1 


.34 


-1.26 



TABLE 7 

Simulations for Series F 



Run 


Mtotal[lO"M0] 


-'Vtotal 


^ merge 


tides 


^ accr 


^halo 
tides 


N ■ * 

^ ^eject 


f surv 


Acs 


'^start 


end 


Fl 


1394.6 


1023 


386 


153 


88 


4 


1 


0.382 


0.138 


-1, 


.38 


-1.28 


F2 


1374.4 


866 


304 


96 


47 


1 


2 


0.480 


0.106 


-1, 


.38 


-1.33 


F3 


1431.1 


906 


277 


183 


196 


1 





0.275 


0.144 


-1, 


.39 


-1.25 


F4 


1163.2 


754 


211 


180 


109 








0.337 


0.126 


-1, 


.33 


-1.27 


F5 


1239.5 


788 


288 


166 


199 








0.171 


0.158 


-1, 


.33 


-1.21 


F6 


1158.8 


527 


141 


60 


55 


1 


1 


0.510 


0.092 


-1, 


.34 


-1.29 


F7 


1300.7 


663 


253 


81 


70 


1 


1 


0.388 


0.113 


-1, 


.37 


-1.28 


F8 


1298.4 


773 


295 


124 


108 


2 





0.316 


0.134 


-1, 


.36 


-1.29 


F9 


1299.5 


652 


226 


93 


61 


2 


1 


0.413 


0.153 


-1, 


.39 


-1.38 


FIO 


1169.4 


589 


182 


83 


72 





1 


0.426 


0.133 


-1, 


.37 


-1.33 


Fll 


1145.8 


741 


269 


148 


76 


1 


1 


0.332 


0.154 


-1, 


.38 


-1.29 


F12 


1227.5 


732 


194 


58 


26 





31 


0.578 


0.067 


-1, 


.39 


-1.36 


F13 


1334.1 


952 


322 


175 


85 


1 





0.388 


0.121 


-1, 


.35 


-1.29 


F14 


1286.0 


844 


272 


180 


107 





1 


0.336 


0.197 


-1, 


.34 


-1.30 


F15 


1470.8 


888 


306 


103 


54 


1 





0.478 


0.108 


-1, 


.36 


-1.31 


F16 


1454.4 


966 


314 


155 


135 


2 


1 


0.372 


0.097 


-1, 


.34 


-1.29 



TABLE 8 

Simulations for Series G 



Run 


Mtotal[lO"M0] 


.'Vtotal 


^ merge 


tides 


^ accr 


Arhalo 
"tides 


.'Veject 


/surv 


Acs 


'^start 


"end 


Gl 


721.9 


372 


72 


55 


6 


47 


1 


0.513 


0.486 


-1.27 


-1.26 


G2 


1034.3 


637 


129 


80 


4 


63 


7 


0.556 


0.405 


-1.29 


-1.26 


G3 


821.1 


530 


127 


89 


5 


73 


1 


0.443 


0.587 


-1.31 


-1.29 


G4 


992.3 


457 


113 


61 


9 


44 


1 


0.501 


0.211 


-1.28 


-1.26 


G5 


899.7 


618 


94 


56 


7 


64 


8 


0.629 


0.443 


-1.29 


-1.28 


G6 


947.1 


452 


95 


31 


2 


34 


17 


0.604 


0.263 


-1.31 


-1.32 


G7 


865.2 


542 


170 


91 


7 


77 


1 


0.362 


0.543 


-1.31 


-1.25 


G8 


1011.6 


725 


169 


101 


7 


71 


1 


0.519 


0.502 


-1.24 


-1.20 


G9 


1100.6 


726 


214 


144 


17 


103 


1 


0.340 


0.467 


-1.27 


-1.23 


GIO 


1174.4 


619 


190 


76 


2 


57 


6 


0.465 


0.458 


-1.28 


-1.23 



number of galaxies destroyed by the tidal field of the 
NFW cluster halo. This halo tidal disruption was nil 
(in the Virgo-like cluster) to a handful (in the Perseus- 
like cluster) for the /^-model background halo. With the 
NFW profile, the number of halo tides is comparable to 
the galaxy tides, even exceeding the latter in runs G5 
and G6. 

The reason for such a behavior is that the NFW halo 
mass profile rises much more steeply than the /3-model 



mass profile of a Virgo-like cluster up to a distance r ~ 
1.9 Mpc. So the NFW halo is significantly more massive 
(by factors as high as 4 — 5) than the /3-model halo at 
distances r < 1 Mpc, where halo tides are dominant (as 
discussed in i j4.5p . Consequently galaxies near the cluster 
center experience a larger tidal field and are more likely 
to be tidally disrupted. 

This larger number of halo tides alters several results 
in our simulated NFW model cluster as compared to the 
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TABLE 9 
Simulations for Series H 



Run 






^ merge 


I 

tides 


^ aecr 


tides 


.^eject 


f surv 


Acs 


'-'start 


end 


HI 


894.3 


755 


171 


152 


5 


122 





0.404 


0.620 


— 1.33 


— 1.33 


H2 


866.2 


556 


195 


68 


5 


92 


I 


0.351 


0.563 


— 1.28 


— 1.24 


H3 


955.0 


621 


161 


84 


2 


67 


\ 


0.493 


0.396 


— 1.33 


— 1.35 


H4 


758.6 


588 


153 


113 


7 


109 


1 


0.349 


0.702 


-1.30 


-1.30 


H5 


1103.6 


997 


197 


132 


5 


68 


22 


0.575 


0.375 


-1.30 


-1.26 


H6 


1008.7 


774 


266 


162 


10 


148 





0.243 


0.682 


-1.29 


-1.24 


H7 


1027.4 


639 


101 


40 


4 


47 


28 


0.656 


0.316 


-1.32 


-1.30 


H8 


833.8 


601 


126 


75 


5 


70 


7 


0.529 


0.496 


-1.29 


-1.28 


H9 


740.3 


494 


109 


53 


9 


54 


4 


0.536 


0.398 


-1.30 


-1.28 


HIO 


952.9 


744 


159 


96 


13 


92 


1 


0.515 


0.451 


-1.35 


-1.32 


Hll 


779.4 


519 


113 


54 


9 


59 


2 


0.543 


0.375 


-1.31 


-1.25 


H12 


993.6 


741 


163 


86 


3 


57 


10 


0.570 


0.306 


-1.33 


-1.29 


H13 


931.8 


534 


131 


60 


5 


62 


2 


0.513 


0.412 


-1.31 


-1.33 


H14 


1012.8 


748 


202 


117 


8 


95 


1 


0.434 


0.508 


-1.30 


-1.26 



TABLE 10 
Simulations for Series I 



Run 


Aftotal[lO"M0] 


-'Vtotal 


^ merge 


tides 


^ aecr 


^halo 
tides 


N ■ * 

Reject 


/ surv 


Acs 


'^start 


end 


11 


1189.1 


495 


114 


53 


26 


75 


1 


0.457 


0.322 


-1.34 


-1.31 


12 


1255.5 


705 


183 


100 


16 


87 





0.452 


0.326 


-1.34 


-1.32 


13 


1189.9 


643 


152 


127 


35 


119 





0.327 


0.391 


-1.29 


-1.31 


14 


1153.5 


548 


108 


78 


8 


66 


1 


0.524 


0.255 


-1.32 


-1.30 


15 


1309.4 


821 


228 


139 


19 


87 


1 


0.423 


0.341 


-1.33 


-1.31 


16 


1330.8 


596 


105 


56 


2 


54 


9 


0.621 


0.265 


-1.31 


-1.30 


17 


1315.8 


775 


215 


126 


37 


105 





0.377 


0.299 


-1.33 


-1.31 


18 


1188.1 


691 


118 


91 


12 


71 


1 


0.576 


0.248 


-1.36 


-1.35 


19 


1236.6 


574 


103 


73 


12 


56 


1 


0.573 


0.295 


-1.32 


-1.30 


no 


1093.4 


634 


96 


96 


8 


64 


1 


0.582 


0.274 


-1.31 


-1.29 



/3-niodel. The mergers exceed the galaxy tides, usually 
by factors 1.3-1.8 (except runs G6 and GIO, where the 
factors are 3 and 2.5). But when tides by galaxy and 
cluster halo are added together, they become comparable 
to or even exceed the number of mergers. The accretions 
are always small in number, and when added to mergers 
do not have much effect on the above. 
The mean /ics for series G is 

/jcs = 0.436 ±0.117. (35) 

This is significantly larger than the ICS mass fraction 
obtained with the /3-model clusters The reason is again 
the numerous halo tides. Some massive galaxies are being 
destroyed by the tidal field of the NFW halo, when they 
come near the cluster center, and this is contributing a 
large mass fraction to the ICS. 
In this series G, we obtained the average values of a 

as 

a,tart = -1-285 ±0.022, (36) 

a^„d = -1-258 ±0.034. (37) 

Also combining the set of galaxies of all the 10 runs in 
series G, we obtained the best fit Schechter exponent for 
the surviving galaxy distribution at z = as a = —1.25. 
Analogous to our approach for the /3-model in §4.3( we 
note that \a\ decreases by ~ 0.03 between z = 1 and 
z = 0. So we performed a new series of simulations, 
series H, using Ostart = —1.31, expecting that this will 
evolve to a ~ —1.28 at z = 0. This series includes galaxy 
harassment but no cD galaxy. The results for series H 
are shown in Table [9l 



Series H continues the trends of series G pertaining to 
a NFW profile. There are a large number of halo tides 
that dominate the mass fraction, and result in a high 
value of /ics • The combined numbers of tidal disruptions 
(by galaxy and halo) are comparable to or exceed the 
numbers of mergers. The mean /ics for series H is 

/ics = 0.471 ±0.129. (38) 

In this series H, we obtained the average values of a 

as 

a,tart = -1-310 ± 0.020, (39) 
a<,„d = -1-288 ±0.036. (40) 

We then performed a series of simulations by putting 
a cD galaxy at the center of the NFW cluster halo, and 
called it series I. The results are shown in TableHU] Here, 
the numbers of tides by the cluster halo and by other 
galaxies are comparable; when added the total occur- 
rence of tides compares to or exceeds that of mergers. 
Comparing series H and series I (NFW-type clusters re- 
spectively without and with a cD galaxy,) there are more 
accretions when a cD galaxy is introduced (similar to se- 
ries D and F). The trend seen before with the Perseus- 
like clusters (between series E and F), that the number of 
tidal disruptions by the background halo reduces in the 
presence of a cD galaxy, is almost absent in the NFW- 
type clusters. Galaxies approaching the cluster center 
get destroyed by the tidal field of the halo before the cD 
galaxy can have any effect. 

The galactic mass fractions dispersed into the ICM are 
neither too high, nor too low. We suspect this is the 
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TABLE 11 
Series for Parameter Variations 
OF /3-MoDEL Halo Density Profile 



Series 


/3 


Pa [g cm" 




r'c [kpc] 


Bbl 


0.3 


1.0 X IQ- 


26 


50 


Bb2 


0.4 


1.0 X IQ- 


26 


50 


Bb3 


0.5 


1.0 X IQ- 


26 


50 


Bb4 


0.6 


1.0 X IQ- 


26 


50 


Bb5 


0.8 


1.0 X IQ- 


26 


50 


Bb6 


1.0 


1.0 X IQ- 


26 


50 


Brl 


0.5 


1.0 X IQ- 


26 


10 


Br2 


0.5 


1.0 X IQ- 


26 


100 


Br3 


0.5 


1.0 X IQ- 


26 


200 


Br4 


0.5 


1.0 X IQ- 


26 


300 


Br5 


0.5 


1.0 X IQ- 


26 


400 


Br6 


0.5 


1.0 X IQ- 


26 


500 



combined effect of putting a cD galaxy in a NFW type 
cluster. There is a tendency of the ICS mass fraction 
to be high in a NFW model cluster, and a cD galaxy 
tends to reduce the mass fraction imparted to the ICM. 
These two opposing trends cause the /ics values to be 
moderate in series I. Here the mean /ics is 



/j(.g = 0.302 ± 0.044 . 



(41) 



5. DEPENDENCE ON PARAMETERS OF CLUSTER HALO 
DENSITY PROFILE 

5.1. (3-Model 

We investigated the dependence of the galaxy out- 
comes and the ICS mass fraction on the parameters gov- 
erning the cluster halo density profile. For the /3-model 
halo density (see ^2.3p . we found that the typical values 
of the relevant parameters cover the range /3 = 0.3 — 1, 
the core radius Tc = 3 — 600 kpc, the core number density 
no = (0.1 — 100) X lO^'^ cm^'^, corresponding to a physical 
density po = (0.017- 17) x lO'^^ g cm'^. W e obtained 
these va lues from several obser v ationa l studies. iLea et al. 
1973[) : lAbramopoulos fc Kul (fl983D: iJones fc Forman 



1984); 



Makino. 



Waxman 
Sasaki, fc Suto 



Miralda-Escud d" 
(1998); Cir ardi et all 



1995 ): 



19981): 



Ettoril (I200dl): IXue fc Wul (i2000D: jBiyiano fc GiraTdi 
(120031) : iPiffaretti fc Kaastral (|2006l) : iMaughan et al 
(120071 ). 



We performed a series of simulations, as shown in Ta- 
ble [Til doing 5 runs in each series, for a total of 60 sim- 
ulations. These were done using a Schechter mass dis- 
tribution exponent a = —1.36 (see i 34.3p . and including 
galaxy harassment, but not including a cD galaxy. 

In series Bbl - Bb6 (Table [TT|) . we explored values 
of j3 between 0.3 - 1.0, with fixed Tc — 50 kpc, and 
po = 1-0 X 10~^^ g cm~^ (which corresponds to a parti- 
cle number density no = 6.0 x 10"'^ cm~'^). We found the 
arithmetic means and the uncertainties on the arithmetic 
means (= cr/>/5) of the galaxy fraction results from the 
5 runs in each series. The resulting average fractions of 
each outcome are plotted in Figure [5] as a function of 
/3, with the galaxy number fractions in the upper panel, 
and the mass fractions in the lower panel. Occurrences 
of ejection out of the cluster f H2.4.5p are not shown, since 
they are either zero or negligibly small compared to the 
other outcomes. As a note, for a certain value of /3, the 
number fractions of the different galaxy outcomes will 
add up to 1. But the added mass fraction of all the 
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Fig. 5. — Number fractions (upper panel) and mass fractions 
(lower panel) of galaxies having different outco mes a s a function of 
parameter /3 of /3-model halo density profile (sec JSTTJ. The symbols 
show the mean galaxy fractions and the error bars are the errors on 
the mean values. The common plotting symbols for both panels 
are: diamond - merger, square - tidal destruction by a massive 
galaxy with fragments going to ICM, filled circle - accretion after 
galaxy tidal disruption, triangle - tidal destruction by the halo, 
open circle - survival up to the present. In addition, the lower panel 
shows the galactic mass fraction that ends up in the ICM, /ics, 
with asterisk plotting symbol, and joined by a dashed line. On this 
panel, the asterisks and the dashed line are indistinguishable from 
the squares and the line connecting them. 



outcomes will be > 1, because there is a double count- 
ing of the mass of galaxies merging and accreting. The 
mass fraction of galaxies surviving, the mass imparted to 
the ICM by tidal disruptions (due to other galaxies and 
cluster halo), and that ejected, will add up to 1. 

From Figure [5] we see that an increase in (3 causes 
a decline of galaxy interactions, leading to survival of 
more and more galaxies up to ^ = 0. The numbers of 
mergers dominate over galaxy tides by a factor of 2 or 
more. But only the low-mass galaxies merge, making the 
merged mass fraction significantly smaller than the tidal 
destruction mass fraction when j3 < 0.6. With increase 
of P from 0.3 to 0.6, destruction by galaxy tides (with 
the fragments imparted to the ICM) decline dramatically, 
mergers have a slight decrease in number but increase in 
mass (implying that more massive galaxies are merging) , 
accretions increase but always remaining below mergers 
and galaxy tides. In between /? = 0.6 — 1.0 there is not 
much change in the resulting galaxy fractions, the trends 
continue slowly as f3 reaches 1. A noteworthy result is 
that at /3 > 0.6 galaxy accretion dominates the mass 
fraction compared to the other outcomes. Tidal disrup- 
tion by the halo and ejection out of the halo arc cither 
zero, or negligibly small in number and mass, for all /?. 
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200 300 
fc (kpc) 

Fig. 6. — Number fractions (upper panel) and mass fractions 
(lower panel) of galaxies having different outcomes, and mass frac- 
tion imparted to the ICM, /ics (lower panel), as a function of 
the core radius, Tc, of /3-model halo density profile. The plotting 
symbols are same as in Figure [5] See i|5.1l for discussions. 



The mass fraction /icg of galaxies imparted to the ICM 
by tidal disruptions is shown in the lower panel of Fig- 
ure [H] with asterisk symbols joined by a dashed line. This 
curve is not distinguishable since, for this set of simula- 
tions, the ICS mass comes only from the dwarf galaxies 
destroyed by the tides of more massive galaxies. So /ics 
coincides with the solid line showing the galaxy tides 
fraction, with the asterisk symbols coinciding with the 
squares. As /3 rises between 0.3 — 0.6, the ICS mass frac- 
tion decreases from ^ 0.3 to ~ 0.1. After that /ics is 
essentially constant at ~ 0.1 when /? > 0.6. 

Next, we explore variations of the core radius rc in 
between 10 - 500 kpc, in series Brl - Br6 (Table fTTj) . 
Here (3 = 0.5, and po = 1.0 x 10~^^ g cm~^ are kept 
fixed. We calculated and analyzed the results in a way 
analogous to that done for series Bbl - Bb6. The means 
and the errors on the means (from the 5 random runs 
in each series) of the resulting galaxy fractions with a 
specific outcome are plotted in Figure [H] as a function of 

From Figure \6\ the following trends can be inferred 
about the dependence of the results on the core radius. 
As Tc rises, a decreasing number of galaxies survive up to 
z = 0, leading to an increase in the mass incorporated to 
the ICM. When rc > 50 kpc, the mass fraction is increas- 
ingly dominated by galaxies getting destroyed by other 
galaxy tides with the fragments being dispersed into the 
ICM. Mergers outnumber galaxy tides when rc is small, 
become comparable to tides at rc ^ 350 kpc, and at 
large-rc tides dominate. But only the low-mass galaxies 



merge (similar to series Bbl - Bb6), causing the merged 
mass fraction to be significantly smaller than that tidally 
destroyed. As rc increases, galaxy tides increase substan- 
tially by number and mass, mergers remain almost con- 
stant with a small reduction at large-rc, and accretions 
(which have a small contribution) decline. There are a 
few cases of tidal destruction by the cluster halo when 
rc > 100 kpc, which has a few % contribution by number 
and mass. The number of occurrences of ejection from 
the halo is always negligibly small in number and mass. 

The lower panel of Figure El shows /icSj the mass frac- 
tion transferred to the ICM by destruction due to tides, 
as asterisk symbols joined by a dashed line, /ics is largely 
dominated by the galactic mass destroyed by the tides 
of more massive galaxies (squares joined by solid lines). 
There is a small contribution coming from tidal disrup- 
tions by the cluster halo, which gets added to the resul- 
tant /ics- With increase of rc between 10 — 500 kpc, the 
ICS mass fraction rises from ~ 0.05 to ~ 0.45. 

5.2. NEW Model 

We explored the dependence of the galaxy outcomes 
and the ICS mass fraction on the parameters gov- 
erning the NFW model cluster halo density profile 
(see ij2.3[) . Typical ranges of the profile parame- 
ters cover values of the scale radius = 100 — 500 
kpc, and the concentration parameter c = 3 — 6. 
We obta ined these valu es from s everal observational 
works , iCarlberg et al.l (I1997D: Ivan der Marel et"aLl 
|[200 0l): lArabadiis Baut z. fc Garmire' (200^; 



Prat t fc ArnaudI (120021)" iBiviano fc Girardi ( 200^); 
Point ecouteau. Arnaud fc Pratt! |2005f ) ; IPratt fc ArnaudI 
(2005i): lMaughan et all (|2007| ). 

Table [T^ shows the series of simulations we performed, 
doing 5 random runs in each series, for a total of 35 
simulations. For this set, we used a Schechter mass dis- 
tribution exponent a = —1.31 (see reasoning in M.6|l . 
included galaxy harassment, but did not include a cD 
galaxy 

In series Nrl - Nr5 (Table [T2|). we investigate five dif- 
ferent values of the scale radius within 10 - 300 kpc, 
keeping the concentration fixed at c = 4.5. We present 
the results in a manner analogous to that done for series 
Bbl - Bb6 and Brl - Br6 in gSH Figure [7] shows the 
means and the errors on the means (from the 5 random 
runs in each series) of the resulting galaxy fractions with 
a specific outcome as a function of r^. 

Some general trends are clear from Figure [7] about the 
dependence of the results on the scale radius. With in- 
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Fig. 7. — Number fractions (upper panel) and mass fractions 
(lower panel) of galaxies having different outcomes, and mass frac- 
tion imparted to the ICM, /ics (lower panel), as a function of the 
scale radius, r^, of NFW-model halo density profile. The plotting 
symbols are same as in Figure [5] See i|5.2l for discussions. 



crease of r^, a smaller number of galaxies survive up to 
the present, causing greater mass transferred to the ICM. 
The most noteworthy feature of these NFW simulations 
is that a large galactic mass fraction is tidally destroyed 
by the cluster halo. In fact when Vg > 100 kpc, the mass 
fraction is increasingly dominated by galaxies getting de- 
stroyed by tidal field of the halo, with the fraction reach- 
ing as high as 0.6 when = 300 kpc. These disrupted 
galaxy fragments are unconditionally dispersed into the 
ICM, which increases the ICS mass fraction substantially. 

The number of mergers is higher than that of tides 
by galaxy and halo when Vs is small; but these numbers 
of mergers, galaxy tides and halo tides become compa- 
rable at rs ~ 300 kpc. At the same time, the galax- 
ies merging have smaller masses (similar to the results 
in i jS.ip . making the merged mass fraction comparable 
to that tidally destroyed by galaxies, both of which are 
hugely outweighed by halo tides. As increases, tidal 
disruption by the halo rise significantly, galaxy tides in- 
crease slightly in number but decrease in mass (implying 
lower mass galaxies are tidally destroyed by other galax- 
ies), and mergers decline. Similar to the results in i j5.ll 
ejection out of the cluster is always negligibly small in 
number and mass. 

The galactic mass fraction imparted to the ICM by 
tidal destructions is shown in the lower panel of Figure [71 
As discussed in previous paragraphs, here /ics is largely 
dominated by the galactic mass destroyed by the tides 
of the cluster halo {triangles joined by solid lines). Tidal 
disruptions by other galaxies have a small contribution 
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Fig. 8. — Fractional number of galaxies destroyed by mergers, 
/d^sJfoJed (fills'^ circles), and by tides, (open circles), 

averaged over all runs within each series of Table [T] Error bars 
show the standard deviation. 



to /ics- As Ts rises between 10 — 300 kpc, the ICS mass 
fraction grows from ~ 0.1 to ~ 0.65. 

Finally, we performed 2 series of simulations, Ncl and 
Nc2 (Table fT2)) . varying the concentration parameter to 
c — A and 6, with a fixed value of the scale radius 
rs = 100 kpc. We found that the mergers continue to 
outnumber the tides. The average values of the ICS mass 
fractions are, /ics = 0.196 with c = 4, and /ics — 0.278 
with c = 6. 



6.1 



6. DISCUSSION 

Tides 



Mergers vs. 

To quantify the relative importance of destruction by 
tides and by mergers, we calculated, for each run, the 
following fractional numbers: 



/•mergers 

destroyed at 

^destroyed 



/■tides 

J destroyed ' 



^t1des+«= 



tides 



where A^, 



destroyed 



A^dcstroyed 



Arhalo 
tides- 



(42) 

(43) 
We 



then averaged the fractions over all the runs in each se- 
ries of the simulations. The results for the set of series 
in Table [T] are shown in Figure [H The destruction by 
mergers clearly dominates over destruction by tides for 
the f3 model, while they are of comparable importance 
for the NFW model. 

Figure [H gives the destruction fraction results for the 
set of series in Tables [11] and [12] For greater values of 
the index /3 of the /3-model density (series Bbl-Bb6, Ta- 
ble [Tl]), mergers rise and tides decline, causing mergers 
to increasingly dominate over tides. The trend is oppo- 
site for increasing core radius Tc of the /9-model (series 
Brl-Br6, Table [TT|) . and for increasing scale radius of 
the NFW model (series Nrl-Nr5, Table [T^: here mergers 
reduce and tides grow, such that finally (at Tc > 400 kpc, 
and Ts > 300 kpc) tides outnumber mergers. Destruc- 
tion by mergers well dominates that by tides in series 
Ncl and Nc2 (Table [Hj). 

6.2. Intracluster Stars 
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Fig. 9. — Fractional number of galaxies destroyed by mergers, 
^mergers (fjUg^j circles), and by tides, fi"*^'^ , (open circles), 

•> dcstroyca ' ' ^ ' destroyed^ ^ ' ' 

averaged over all runs within each scries of Tables [TTI and [121 Error 
bars show the standard deviation. 



Several mechanisms can contribute to removing stars 
from individual galaxies and putting them in the intra- 
cluster space. The efficiency and relative importance of 
these processes is expected to vary according to the lo- 
cation inside a cluster and during its evolution. If ram- 
pressure stripping or harassment are dominant mecha- 
nisms to produce IC stars, then the ICL fraction should 
increase with the mass of a cluster. On the other hand, if 
galaxy-galaxy merging is the dominant mechanism, and 
most of the ICL formed early on in cluster collapse, then 
the ICL fraction should be independent of present cluster 
mass. 

The ICL fraction depends on the merger history of 
the cluster, presence of a cD galaxy, and the mor- 
phology of the cluster galaxies. Observations (e.g., 
iKrick fc Bernsteinl 120071 ) show that if a cD galaxy is 
present then the ICL profile is centrally concentrated, im- 
plying that the ICL formed by galaxy interactions at the 
center, or formed earlier in protoclusters and later com- 
bined at the center. The ICL fraction should evolve with 
redshift, as the number of galaxy interactions increase 
with time. Cosmological simulations indicate that the 
ICL fraction does in c rease with time as clusters evolve 
(jWillman et al.l[200l iR^dick et al.l[2006l ). 

The different formation mechanisms impart some dis- 
tinct perceptible trends in the ICL. If most of the IC stars 
originate in initial cluster collapse, their distribution and 
kinematics should closely follow that of the galaxies in 
the cluster. If the ICL builds up slowly with time due 
to processes like "galaxy harassment," "tidal stripping," 
etc., then a fraction of IC stars should be located in long 
streams along the orbits of parent galaxies. Cosmolog- 
ical simulations with very high resolution is required to 
address these issues numerically. 

There have been several observational measurements 
of the light (or, mass) fraction contained in the ICS with 
respect to the total light in a cluster. We collected some 
values of the ICS fraction from the literature, and list 
them in Table [T31 In Figure [TO] we show the ICS mass 
fraction we obtained in our simulations, plotted as hor- 
izontal lines showing the average /ics from the runs in 
Series A-I (Table [1] For comparison, the observed 
/ics values (from Table I13p are shown by the symbols 
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Fig. 10. — Fraction /jQg of intracluster stars. The horizontal 
lines show the average /ics values in our simulations, from the 
runs of the series in Table [T](j|4]l, with solid line: Virgo-like cluster 
(series A— D), dashed line: Perseus-like cluster (series E-F), dotted 
line: NFW model cluster (series G-I). The symbo ls a nd error bars 
show actual measurements, as tabulated in Table [131 



and error bars. We can clearly see that the ICS mass 
fraction in clusters from observations fall well within our 
simulation predictions. A few clusters have too small 
/ics 7 which are probably galaxy groups and low-mass 
clusters. 

Our simulation results indicate that the tidal destruc- 
tion of dwarf galaxies (by other galaxies and by the 
cluster halo) in clusters can sufficiently explain the ob- 
served fraction of intracluster light. Supporting our 
conclusions, observational studies do find a significant 
ICL component arising from dwarf galaxies in clusters. 
From spectroscopic and photometric analvsis. lCote et"aLl 
(|2001h found that the metal poor globular cluster sys- 
tem associated with M87 (the cD galaxy in Virgo) is 
not formed i n situ, but is s t rippe d from the dwarf galax- 
ies in Virgo. iDurrell et al.l ()2002f ) observed that most of 
the IC red-giant branch stars in Virgo are moderately 
metal rich suggesting that they are tidally stripped from 
intermediate-luminosi ty galaxies, thoug h strip ping from 
DCs is not ruled out. iMcNamara et al.l (|1994D observed 
trails of of IC stars and gas from the tidal disruption 
of a DG by a companion giant elliptical galaxy in Virgo 
cluster. 

Results from our simulation runs (in 21 SJS]) indi- 
cate that, for each cluster halo density profile, namely, 
f3, and NFW models, /ics increases with the mass of the 
cluster halo. This is consistent with studies finding that 
more massive clusters have a larger fraction o f ICL than 
the le ss massive ones |Lin fc Mohr.,2004 : .Murante et alj 
[200l . 

6.3. Limitations of the Method 

The strengths and weaknesses of the methodology used 
in this work both reside in our somehow original ap- 
proach of using one single particle to represent each 
galaxy, combined with a subgrid treatment of galaxy 
mergers, tidal disruption, and galaxy harassment. 

On the positive side, this approach has enabled us to 
perform a very large number of simulations (227 total), 
covering a fairly large parameter space, while obtaining 
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TABLE 13 

Observed values of the ICL mass as a fraction of the 
total cluster mass 



Index 


Cluster 


/ICS (%) 


A/ics 


1. 


Coma 


50 




2. 


Abell 1689 


30 




3. 


Abell 1651 


< 5 




4. 


M96 (Leo) Group 


< 1.6 




5. 


HCG 90 


45 


5 


6. 


Virgo 


15.8 


8 


7. 


A801 


16 


4.7 


8. 


A1234 


17 


4.4 


9. 


A1553 


21 


16 


10. 


A1914 


28 


16 


11. 


93 clusters 


50 


10 


12. 


683 clusters 


10.9 


5.0 


13. 


A4059 


22 


12 


14. 


A3880 


14 


6 


15. 


A2734 


19 


6 


16. 


A2556 


6 


5 


17. 


A4010 


21 


8 


18. 


A3888 


13 


5 


19. 


A3984 


10 


6 


20. 


A141 


10 


4 


21. 


AC 114 


11 


2 


22. 


AC 118 


14 


5 



Reference 



Bernstein et al. (1995) 
Tyson & Fischer ( 1995) 
^onzalez et al. (2000J 
Castro-RodriKuez et al. (2003) 
White et al. (2003) 
FV;ldmeier et al. (2004b) 
Rjldmeier et al. ( 2004a) 
F\3ldmeier et al. ( 2004a) 
Fbldmeier et al. ( 2004a) 
Rildmeier et al. (2004a) 
Lin & Mohr (2004) 
Zibetti et al. (2005) 
Krick & Bernstein (2007) 
Krick & Bernstein (2007) 
_Krick & Bernstein (200T) 
_Krick & Bernstein (200T) 
Krick & Bernstein (20(S 
Krick & Bernstein (20OT) 
Krick & Bernstein (2007) 
_Krick & Bernstein (2007) 
Krick & Bernstein (2007) 
Krick fc Bernstein (2007) 



^ IKrick &c Bernstei measured these /ics values in the r band. 



statistically significant results. In particular, the disper- 
sion of the results seen in Tables B] [TUl i ust ifies a posteriori 
the number of simulations we have performed. Doing this 
many simulations without resorting to subgrid physics 
would have been computationally prohibitive. In addi- 
tion, this approach greatly facilitates the analysis and 
interpretation of the results. In a simulation in which a 
galaxy would be represented by a large numbers of par- 
ticles, a dwarf galaxy could be partly destroyed by tides, 
with some fragments being dispersed in the ICM, other 
fragments accreting onto the larger galaxy, with enough 
of the galaxy surviving that it can still be counted as a 
galaxy. It would then become more difhcult to quantify 
objectively the relative importance of mergers and tidal 
disruptions. 

In implementing the subgrid physics, we have at- 
tempted to make the most reasonable choices possible. 
One free parameter is the geometric factor entering in 
equation (|10p . but for reasonable density profiles, the 
values of that factor do not appear to vary much. The as- 
sumption that a galaxy is considered "tidally disrupted" 
if 50% of its mass becomes unbound is also the most 
reasonable one we could make. 

Our technique for generating the initial conditions is 
based on four key assumptions: (1) the galaxy distribu- 
tions are isotropic, (2) the galaxy number density pro- 
file v{r) follows the density profile phsdoir) of the back- 
ground cluster halo, (3) the mass is segregated in the 
cluster, with the most massive galaxies being located in 
the center, and (4) the cluster is in equilibrium. So even 
though our prescription for generating the initial condi- 
tions contains many tunable parameters, we believe that 
the underlying approach is sound. 

On the negative side, two particular aspects of the 
methodology can be considered weak. First, the treat- 
ment of galaxy harassment is highly speculative. We 
have assumed that some amount of orbital kinetic en- 
ergy AE is dissipated into internal energy during an en- 



counter between two galaxies, that this amount is related 
to the initial internal energies of the galaxies, and that 
the energy dissipated is distributed equally between the 
two galaxies. The dissipation of energy and its conse- 
quences during a real galactic encounter are certainly 
much more complex, so the subgrid model could poten- 
tially be refined. But this would require an intensive 
study of high-speed encounters, with detailed numerical 
simulations, and such study remains to be done. 

Another important limitation of our approach is that 
it deals with isolated clusters in equilibrium. In the real 
universe, clusters constantly experience mergers and ac- 
cretion. We justify our approach by the fact that most 
clusters will experience, at some epoch, a major merger^ 
during which most of the final mass of the cluster is as- 
sembled. From that point, if we can neglect the addition 
of mass by minor mergers and accretion, the cluster can 
be treated as isolated. Of course, such scen ario cannot 
descr ibe all clusters. In a forthcoming paper (jBrito et al.l 
I2008D . we will present a study of cluster formation and 
evolution inside a cosmological volume containing many 
clusters. This will be achieved by implementing the sub- 
grid approach described in this paper into a cosmological 
N-body algorithm. 

7. SUMMARY AND CONCLUSION 

We have designed a simple model for the evolution of 
galaxies in an isolated cluster, in order to compare the 
destruction of dwarf galaxies by mergers vs. tidal disrup- 
tions, and to predict the contribution of dwarf galaxies 
to the origin of intracluster stars. Our algorithm com- 
bines a direct N-body computation of gravitational inter- 
actions, along with a subgrid treatment of the other phys- 
ical processes (merger, tidal disruption, accretion etc.) 
of the galaxies. Using this algorithm, we have performed 
a total of 227 numerical simulations of galaxy clusters, 
examining the fate of dwarf galaxies. Our results and 
conclusions are as follows. 
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(1) The destruction of dwarf galaxies by mergers domi- 
nates over destruction by tides, in most of our simulation 
runs with all the models (/3- Virgo, /9-Perseus, NFW) of 
cluster halo density. The two destruction mechanisms 
become comparable, and tides outnumber mergers, for 
the NFW and the f3- models when their scale/core ra- 
dius approach and exceed ^ 200-300 kpc. 

(2) The destruction of dwarf galaxies by the tidal field 
of other galaxies and by the cluster halo imparts a sig- 
nificant amount of galactic mass into the ICM. This is 
sufficient to account for the observed fraction of intra- 
cluster light in galaxy clusters. In our simulations, the 
ICS mass fraction, /ics, has a range 4.5% — 66.7% (with 
the majority of values being in between 10% — 45%), for 
the different parameter sets of the (3 and NFW models of 
cluster halo density we considered. We see a clear trend 
of increase of /ics with the mass of the cluster halo. All 
these are well consistent with observations and other nu- 
merical studies. 

(3) In the NFW model simulated clusters, there are a 
large number of tidal disruptions of galaxies caused by 
the gravitational potential of the cluster halo, and this 
component dominates the mass fraction. We note that it 



has been our assumption that the cluster halo is station- 
ary, and does not evolve in response to the forces exerted 
on it by the galaxies f ji2.ip . Such an assumption is prob- 
ably a poor one with the NFW model clusters. We point 
out that this could imply a possible solution to the cusp 
crisis of cluster dark matter halos. The central cuspy re- 
gion of the cluster dark matter halo could have inelastic 
encounters with the member galaxies, which could inject 
energy into the halo and erase the cusp. 

(4) In our simulations, the presence of a cD galaxy 
increases occurrences of accretion, decreases tidal dis- 
ruptions by the cluster halo, and reduces the ICS mass 
fraction. This is opposite to the trend seen from obser- 
vations that /ics is higher in the presence of a cD. 



This work benefited from stimulating discussions with 
L. Edwards and C. Robert. We thank John Kormendy 
for useful correspondence. All calculations were per- 
formed at the Laboratoire d'astrophysique numerique, 
Universite Laval. We thank the Canada Research Chair 
program and NSERC for support. 



APPENDIX 
THE INTERNAL ENERGY OF GALAXIES 

Since we represent galaxies as individual particles, we cannot directly compute their internal energy. We therefore 
need an estimate that can then be used in equation (llOp . We write the potential energy of a galaxy of mass M and 
radius R as 

»' = -^. (Al) 

where ^ is the geometric factor, which depends on the shape and density distribution of the galaxy. For a uniform- 
density sphere, C = 3/5. In our simplified model, we treat galaxies as spheres, but we should certainly not assume a 
uniform density. Instead, any galaxy will be centrally concentrated. The value of ^ will depend on the assumed density 
profile, but we do not expect that dependence to be very strong if we stick with reasonable profiles. So we consider 
the simplest case of an isothermal sphere with a cutoff radius R. The density and mass inside r are given by 



where M 



AirRr^ ' 
Mr 

m{R) is the total mass. The gravitational field is given by 



p{r)- 
m{r) - 



M 



g 



Gm{r) ^ 



GM 

7 

rR 



We integrate this expression, with the boundary condition 4){R) = —GM/R, to get the gravitational potential, 

, GM 



R 



In- 



R 



1 



The potential energy is given by 



W 



(j){r)p{r)(Pr 



GM^ 
~2R^ 



In- 



R 



1 Mr 



GM^ 
R 



(A2) 
(A3) 

(A4) 

(A5) 

(A6) 



Hence, C = 1 for isothermal sphere. Interestingly, this is not much different from the value of 3/5 for a uniform 
sphere. This supports our claim that the sensitivity of C on the density profile is weak. For the kinetic energy, we 
assume that the galaxies are virialized. Hence, K = —W/2, and therefore the internal energy is given by 



U = K + W = — 



GM^ 
'2R 



(AT) 
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